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Streszczenie Jednym z zastosowah metod sredniego pola w fizyce jadrowej jest obecnie 
badanie egzotycznych symetrii jader atomowych. Wiaze si§ to w szczegolnosci z analiza rotacji 
jader wokol osi pochylonej wzgl§dem osi glownych rozkladu masy w modelu Tilted- Axis Crank- 
ing (TAC). Niniejsza praca przedstawia jedne z pierwszych obliczeri TAC wykonanych metodami 
w pehii samozgodnymi. Zastosowano w niej metode Hartree'ego-Focka z dwucialowym oddzia- 
lywaniem efektywnym Skyrme'a. Stworzono program komputerowy pozwalajacy na iamanie 
wszystkich symetrii przestrzennych rozwiazania. Jako pierwsze zastosowanie przeprowadzono 
obliczenia dla pasm magnetycznych w 142 Gd oraz chiralnych w 130 Cs, 132 La, 134 Pr i 136 Pm. 
Ich wyst§powanie zwiazane jest odpowiednio z nowym mechanizmem iamania symetrii sfe- 
rycznej oraz ze spontanicznym lamaniem symetrii chiralnej. Samozgodne rozwiazania w 142 Gd 
potwierdzaja istotna role, mechanizmu shears w tworzeniu calkowitego momentu p§du. Zgodnosc 
z danymi doswiadczalnymi nie jest jednak zadowalajaca, prawdopodobnie wskutek nieuwzgl§d- 
nienia korelacji par lub przeszacowania deformacji. Wyniki w 132 La stanowia pierwszy w pelni 
samozgodny dowod, ze rotacja jadrowa moze przybierac charakter chiralny. Wykazano, ze ro- 
tacja chiralna moze wyst§powac tylko powyzej pewnej krytycznej cz§stosci obrotu. Sprawdzo- 
no tez, ze czlony pola sredniego Skyrme'a nieparzyste wzgl§dem odwrocenia czasu nie maja 
jakosciowego wplywu na wyniki. 

Resume Une des applications recentes des methodes de champ moyen en physique nucleaire 
est l'etude des symetries exotiques du noyau. Cette problematique est reliee, en particulier, a 
l'analyse de la rotation nucleaire autour d'un axe incline par rapport aux axes principaux de la 
distribution de masse dans le modele dit de Tilted- Axis Cranking (TAC). Cette these presente 
l'un des premiers calculs TAC effectues dans le cadre de methodes entierement auto-coherentes. 
La methode Hartree-Fock avec l'interaction effective a deux corps de Skyrme a ete utilisee. Un 
code numerique a ete ecrit qui permet de briser toutes les symetries spatiales des solutions. 
Comme premiere application, des calculs pour les bandes magnetiques dans 142 Gd et pour les 
bandes chirales dans 130 Cs, 132 La, 134 Pr et 136 Pm ont ete effectues. L'existence de ces bandes 
est due a un nouveau mecanisme de brisure de la symetrie spherique, et de brisure spontanee de 
la symetrie chirale, respectivement. Les solutions auto-coherentes dans 142 Gd confirment le role 
important du mecanisme shears dans la generation du moment angulaire. Pourtant, l'accord 
avec les donnees experiment ales n'est pas satisfaisant, probablement a cause de l'omission des 
correlations d'appariement dans les calculs ou de la possible surestimation de la deformation. 
Les resultats obtenus dans 132 La constituent la premiere preuve entierement auto-coherente que 
la rotation nucleaire peut acquerir un caractere chiral. II a ete demontre que la rotation chirale 
ne peut avoir lieu qu'au-dessus d'une certaine frequence angulaire critique. II a ete egalement 
verifie que les termes du champ moyen de Skyrme impair par rapport au renversement du temps 
n'ont pas d'influence qualitative sur les resultats. 

Abstract Currently, one application of the mean-field methods in nuclear physics is the 
investigation of exotic nuclear symmetries. This is related, in particular, to the study of nuclear 
rotation about an axis tilted with respect to the principal axes of the mass distribution in the 
Tilted- Axis Cranking (TAC) model. The present work presents one of the first TAC calculations 
performed within fully self-consistent methods. The Hartree-Fock method with the Skyrme 
effective two-body interaction has been used. A computer code has been developed that allows 
for the breaking of all spatial symmetries of the solution. As a first application, calculations 
for the magnetic bands in 142 Gd and for the chiral bands in 130 Cs, 132 La, 134 Pr, and 136 Pm 
have been carried out. The appearance of those bands is due to a new mechanism of breaking 
the spherical symmetry and to the spontaneous breaking of the chiral symmetry, respectively. 
The self-consistent solutions for 142 Gd confirm the important role of the shears mechanism in 
generating the total angular momentum. However, the agreement with experimental data is 
not satisfactory, probably due to the lack of the pairing correlations in the calculations or to 
the possibly overestimated deformation. The results obtained for 132 La constitute the first fully 
self-consistent proof that the nuclear rotation can attain a chiral character. It has been shown 
that the chiral rotation can only exist above a certain critical angular frequency. It has also been 
checked that the terms of the Skyrme mean field odd under the time reversal have no qualitative 
influence on the results. 
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Chapter 1 
Introduction 



The phenomenon of spontaneous symmetry breaking is a leading theme of a multitude of 
quantum effects. For instance, it is due to the spontaneous breaking of the particle-number 
symmetry that superconducting condensates appear in metals. Spontaneous breaking oc- 
curs if a system in its endeavor to attain the minimal energy chooses a symmetry-violating 
state even though the underlying interactions are invariant under the concerned symme- 
try. Nevertheless, it is the nature of the interactions that determines which symmetries 
are broken and under which conditions. Therefore, study of symmetry-violating states 
brings one closer to understanding the interactions - the most fundamental goal in physics. 

In the nuclear structure physics, one important example of spontaneous symmetry 
breaking is the existence of deformed nuclei. Like for molecules, violation of the spherical 
invariance leads to the appearance of rotational excitations, which manifest themselves 
in specific sequences of levels, the rotational bands. Depending on the conservation or 
violation of other symmetries, like the plane reflection, those bands can have different 
structures. Recently, two novel types of bands have been observed experimentally, that 
have challenged the theoretical models. They are called magnetic and chiral bands. 

The peculiarity of the magnetic structures is that, contrary to all the previously known 
rotational bands, they appear in nuclei whose charge distribution is nearly spherical, which 
can be inferred from the very weak electric quadrupole transitions. On the other hand, 
magnetic dipole transitions, to which the bands owe their name, are strong. Thus, the 
bands constitute the first evidence for the breaking of the spherical symmetry by a large 
dipole moment, which, in turn, is produced by highly asymmetric current distributions. 
In the absence of charge deformation, no collective rotation is possible, and the magnetic 
bands also entail a new mechanism of generating the angular momentum, in which valence 
nucleons play a crucial role. The angular momenta of the valence particles align coherently 
along one direction, and those of the valence holes along a perpendicular direction. The 
spin of the levels in the band is produced by gradual alignment of the two angular momenta 
vectors. This has been dubbed shears mechanism, because it resembles the closing of a 
pair of shears used for cutting the sheep wool. The magnetic bands were first observed in 
early 1990s, and more that 130 such structures have been found so far. 

The chiral bands have the form of doublets of bands lying very close in energy. Since 
other hypotheses failed to reproduce the vanishing energy splitting between them, it has 
been suggested that the doublets are due to the possible existence of two enantiomeric 
(left- or right-handed) forms of the nucleus. The chiral bands are observed mainly in well- 
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deformed nuclei, in which there is one valence proton particle and one valence neutron 
hole in an orbital of high angular momentum. The former drives the nucleus towards 
elongated shapes, while the latter towards oblate ones. The interplay of these opposite 
tendencies may result in a shape resembling a triaxial ellipsoid. In the triaxially deformed 
nucleus, the particle and hole align their angular momenta along the short and long axes 
of the density distribution, respectively. Moreover, the moment of inertia with respect 
to the medium axis is the largest, which favors the collective rotation around that axis. 
Thus, the total angular momentum vector has non-zero components on all the three axes, 
and those component vectors can form either a left-handed or a right-handed system. 
This means, according to the original Kelvin's definition, that the system is chiral, only 
that its chirality is related to the handedness of angular momenta rather than of position 
vectors. First chiral bands were identified about the year of 2000, and ten odd occurrences 
are known now. 

The chiral rotation has been successfully examined within the phenomenological Par- 
ticle-Rotor Model, in which the nucleus is represented by a rigid core with a few valence 
nucleons coupled to it. However, the very concept of both the shears mechanism and 
rotational chirality came from considerations within the mean-field cranking approach, 
which is one of the most fundamental methods in the nuclear structure physics. In this 
approximation, each nucleon is moving independently in a rotating potential that repre- 
sents an averaged interaction with other nucleons. The main point about the mean-field 
description of the two considered phenomena is that it must allow for rotation about 
an axis that does not coincide with any principal axis of the mass distribution; such a 
variant of the cranking model is called Tilted- Axis Cranking (TAC). Before the interest 
in magnetic and chiral bands arose, only rotations about principal axes were considered, 
so addressing the new effects required an upgrade of the existing numerical tools. A 
first TAC code was written by Frauendorf [Fra93a ; it used a phenomenological mean 
potential to describe the properties of the valence particles, and the nuclear liquid-drop 
model to account for some bulk properties of the nucleus. A great deal of experimental 
data, mainly on magnetic bands, has been analyzed by using that code, and generally a 
good agreement was obtained. As far as more fundamental methods are concerned, only 
one magnetic band in 84 Rb has been studied by Madokoro et al. [MadOOa within the 
Relativistic-Mean-Field method. 

Although the phenomenological model of Frauendorf has led to a remarkable suc- 
cess, a more fundamental description requires self-consistent methods, in which the mean 
potential is indeed generated from averaging an effective two-body interaction with the 
nucleonic density. First of all, such an approach is capable to provide a strong test of the 
stability of the proposed shears and chiral configurations with respect to the core degrees 
of freedom. It is also necessary to take into account several important effects, like all kinds 
of polarization of the core by the valence particles and full minimization of the underly- 
ing energies with respect to all deformation variables, including the deformation of the 
current and spin distributions. Application of self-consistent methods to the description 
of the magnetic and chiral rotation is the subject of this PhD dissertation. 

The author developed a code that can perform TAC calculations within the full 
Hartree-Fock (HF) method with the Skyrme effective interaction [RinOOa . It is one 
of the first existing tools of this kind, and allows for a study of many other effects; for 
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instance, it has already been used to calculate the time-reversal violating Schiff moment 
of 225 Ra |Eng03a| in connection with the search for fundamental breaking of the T or CP 
symmetry beyond the Standard Model. In the present work, the first Skyrme-HF TAC 
solutions were found. For the study of the shears bands, the nucleus of 142 Gd was cho- 
sen. The HF results corroborated that an important portion of the angular momentum 
is generated by the shears mechanism, although the shape deformation is non- negligible, 
and the collective rotation of the core is also present. The chiral solutions were looked 
for in four N = 75 isotones, 130 Cs, 132 La, 134 Pr, 136 Pm, and found in the second one. It 
was established that the chiral rotation appears only above a certain critical value of the 
angular frequency. The origin of the critical frequency was explained in terms of a very 
simple classical model, that provides a surprising agreement with the HF results. 

In Chapter El an introduction to the physics of the magnetic and chiral bands is given, 
together with a review of the literature. Chapter El describes the theoretical tool used 
in this work - the Skyrme-HF TAC method. Particular emphasis is put on the issue of 
spontaneous symmetry breaking in rotational bands and on technical differences between 
the self-consistent and phenomenological approaches. One aspect that can be studied 
uniquely within the self-consistent approach is the influence of the HF time-odd densities 
and fields on the shears and chiral solutions. The way this is investigated in the present 
work is also described in that Chapter. Finally, the code hfodd is presented. The 
rotational behavior of the nuclei here under study is determined, to a big extent, by the 
alignment properties of a few valence hn/2 nucleons. In Chapter EJ their basic features 
are established from standard (non- TAC) cranking calculations and from pure symmetry 
considerations. The Skyrme-HF results for the shears band in 142 Gd are presented in 
Chapter Chapter Ogives the planar and chiral solutions for the N = 75 isotones. There, 
the classical model is formulated, the expression for the critical frequency is derived, and 
its values are analyzed. Since the present calculations constitute the first application of the 
Skyrme-HF TAC method, several technical details on how the solutions were obtained are 
also given in that Chapter. Chapter [7| summarizes the main conclusions from the present 
work, and outlines prospects for further research. 
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Chapter 2 

Magnetic and chiral rotations 



A classical rigid body can rotate uniformly only about the principal axes of its inertia 
tensor. However, already Riemann showed .Rie60aJ that in liquids rotation about a tilted 
axis may take place if there is an intrinsic vortical motion. In a drop of nuclear matter, 
the nucleus, vorticity may arise from the presence of unpaired nucleons near the Fermi 
level. The swirling vortices produce the single-particle (s.p.) angular momenta that take 
various directions and add to the collective rotation drawing the total spin away from the 
principal axes. It is the rotation about a tilted axis that breaks the signature and chiral 
symmetries, and gives rise to the magnetic and chiral bands. The two phenomena show 
different ways of coupling the s.p. spins to the collective rotation. 

2.1 Magnetic rotation 

In early 1990's, very regular rotational bands were observed in several light lead isotopes 
Bal92a| ICla92a[ IKuh92aj , and were initially mistaken for superdeformed structures. But 
detailed measurements showed that they had an unusual feature: weak E2 and strong Ml 
transitions. To date, many more similar structures have been identified. They are called 
magnetic dipole or shears bands. This Section recapitulates their main characteristics and 
presents the mechanisms that are believed to lie at their origin. 

Taking into account the vast experimental evidence, one can enumerate the following 
properties that can be considered as a definition of the magnetic dipole bands: 

• Levels in the band form the J + , (J + l) + , (J + 2) + or I~, + (J + 2)~ spin-parity 
sequence. 

• Bands never start at spins less than about 10h, except for the lightest nuclides, and 
approximately follow the regular rotational dependence of E ~ /(/ + 1). 

• The E2 transitions within the band are very weak, with reduced probabilities B(E2) 
typically not exceeding ~ 0.1 e 2 b 2 . In contrast, the values of B(M1) are exception- 
ally large, in the range of ~ 2 — 10/ijv. They exhibit a characteristic fall with 
increasing spin. 

• The ratio jW/B(E2) assume the values of ~ 100 ^ 2 MeV _1 e- 2 b~ 2 that is roughly 
an order of magnitude larger than for normal or superdeformed bands. 
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Figure 2.1: Experimental evidence for the magnetic bands throughout the nuclear chart. 
Nuclei, in which magnetic bands have been found [AmiOlaJ (May 2001), are marked with 
black squares, and the gray background represents the known isotopes Ant 02a (July 
2002). Principal high-j orbitals available for the valence particles and holes are marked. 



So far, more than 130 bands possessing such properties have been observed in more that 
60 nuclides that are marked in Fig. 12.11 They group in four islands situated not far from 
magic numbers. 

If a nucleus exhibits rotational excitations, its spherical symmetry must be sponta- 
neously broken. Before 1990, there was only evidence for breaking the rotational invari- 
ance by a deformation of the charge distribution. But in the shears bands, the very weak 
E2 transitions unequivocally point to an almost spherical shape. Therefore, the main 
question about the newly found bands is what violates the spherical symmetry, if not 
the shape. On the other hand, sequences of levels linked by relatively strong Ml and 
weak E2 transitions appear quite often in nearly spherical nuclei. They originate from 
recouplings of the angular momenta of valence nucleons, and are well understood within 
the spherical shell model. However, these sequences constitute complicated multiplets 
with energy changing abruptly with spin. From this point of view, the peculiarity of the 
magnetic bands consists in their surprisingly regular level spacing. 

Origins of that regularity were sought in a series of shell-model studies. Frauendorf et 
al. investigated the bands observed in light Pb isotopes |Fra96a and in odd In isotopes 
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Figure 2.2: The shears mechanism; from the bandhead a) to the maximum alignment c). 
The angular momentum, j, is generated by gradual alignment of spins of valence particles 
and holes, j p and j h . The corresponding contributions, \i v and fi h , to the magnetic dipole 
moment, fi, are drawn schematically for the case in which the valence particles are protons, 
and the valence holes are neutrons. 

|Fra97bj using the surface-delta interaction. They found that a substantial portion of 
the angular momentum originates from coupling of high-j valence particles in one kind 
of nucleons with high-j valence holes in the other kind. However, regular structures 
were obtained only if low-j orbitals were included in the configuration space. These or- 
bitals provide a certain polarizability of the core, which introduces the long-range, mainly 
quadrupole correlations. These correlations smooth out the E(I) dependence, because 
they involve many nucleons and average their influence. Therefore, the magnetic rotation 
is a collective phenomenon, despite that mainly the valence high-j nucleons contribute 
to the angular momentum 1 . The core polarization plays a crucial role, even though it 
is weak as opposed to that in well deformed nuclei. These features are supported by 
several experimental observations of transitions from irregular to regular structures, as 
the Z and N numbers get away from magic shells, increasing the polarizability. Frauen- 
dorf [Fra97bJ found such a transition in 49^58,60,62,64 isotopes with increasing N, in good 
correlation with the shell-model results. Similar behavior was reported for isotopes of Cd 
and Sn [FraOlaj . In the region of iV = 80, it was found in 64Gd78 i 79 i go |Rza01allLieT)2a1 . 
Inspecting the table of magnetic bands AmiOla , one easily notices that the effect is also 
well pronounced when passing from the isotopes of soHg to the isotopes of 82Pb. In fact, 
the magnetic bands have been observed in close, but never immediate vicinity of doubly 
magic nuclei, see Fig. 12.11 

If the low-j orbitals are excluded from the configuration space in the discussed shell- 
model study, the valence high-j particles and holes rearrange their angular momenta in 
a complicated way from one level to another. The appearance of quadrupole correlations 
with the inclusion of those orbitals causes the valence high-j particles to orient their an- 
gular momenta coherently in the same direction and the holes to behave similarly. Such a 
stretched alignment of several high-j contributions results in long particle and hole angular 
momenta vectors, j p and j h , which are called blades. In general, both kinds of nucleons 



however, to avoid confusion, in the following the term collective rotation is reserved for the standard 
rotation of the deformed mass distribution. 
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can contribute to each blade, but in most cases the particles are protons and the holes 
are neutrons. At the bandhead, the two blades turn out to form an angle of nearly 90° 
(Fig. 12.2b ). Since the proton and neutron effective ^-factors, g n and g u , have opposite 
signs, the perpendicular coupling also gives rise to a large magnetic dipole moment, /x, 
with a significant component, fi±, perpendicular to the total spin, j, (Fig. 12.2b ). This 
is the source of the strong Ml transitions, because the reduced probability B (Ml) is 
proportional to The perpendicular coupling of j p and j has been confirmed by an 
experimental determination of the effective ^-factor in the bandhead of a shears band in 
193 Pb |Chm97a] . In the band, j p and j h align towards each other, retaining their constant 
lengths. In this way, the total spin increases, while decreases, causing the observed 
drop-off in the Ml strength (Fig. 12.2b ). Assuming the geometry shown in Fig. 12.21 Mac- 
chiavelli and Clark Mac98aJ calculated the B(M1) values for the magnetic bands in 198 Pb 
and 199 Pb as proportional to and obtained a good overall agreement with experiment. 
Obviously, the maximum spin that can be generated in such a way approximately equals 
the sum of j p and j h , whereas [i± and B (Ml) vanish when the maximum alignment is 
achieved (Fig. 12.2b ). This mode of generating the angular momentum and the magnetic 
moment was originally proposed by Frauendorf |Fra93al . and was then dubbed shears 
mechanism [Bal94aJ, because it resembles the closing of a pair of shears used for cutting 
the sheep wool. 

In view of the important role of the quadrupole polarization, the magnetic bands have 
been investigated within approaches that take that polarization into account in a model 
way. In |Mac98bj and als o in pac98allMac98cllMac99aj . Macchiavelli and Clark used the 
particle- vibration theory |Boh75aj . considering the coupling of the valence high-j particles 
and holes to the quadrupole vibrations of the nuclear shape about a spherical equilibrium. 
They postulated the stretched coupling of the constituent angular momenta in each blade 
and treated all the valence particles as one particle with spin j p and all the valence holes as 
one hole with spin j h . In the second order of the perturbation series, the particle-vibration 
coupling gives rise to an effective interaction between the valence nucleons. It takes the 
form of the second-order Legendre polynomial, P2(cosi9), where d is the angle between 
the position vectors of the valence particle and hole. It was shown in |Boh53aj that in the 
limit of high angular momenta the expectation value of the interaction P2(cos , &) in the 
state \j p j h j) becomes proportional to P2(cos9), where, this time, 9 is the angle between 
j p and j h , the shears angle. Actually, these results were obtained in such a regime of the 
particle-vibration theory, in which there already exist a static quadrupole polarization 
of the core. Therefore, it is not surprising that they could also be derived within the 
mean-field approximation. Frauendorf considered |Fra01aj the quadrupole-quadrupole 
model with some simplifications similar to that made by Macchiavelli and Clark. He 
presumed that in both protons and neutrons the ensemble of the valence nucleons can 
be treated as one object with constant quadrupole tensor and mean angular momentum 
pointing along one of its principal axes. Under these assumptions and with the aid of the 
addition theorem for spherical harmonics, the quadrupole-quadrupole interaction between 
the valence particles and holes transforms directly into the form of P2(cos8). He also 
postulated that the valence nucleons are the only source of polarization of the otherwise 
spherical core and that the core energy depends quadratically on the induced quadrupole 
moment. The core defined in this way only modifies the strength of the P2(cos8) force. 
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It comes both from the particle- vibration theory and from the mean-field quadrupole- 
quadrupole model that for the relevant particle-hole case the coupling constant of the 
P2(cos#) interaction is positive. Since the polynomial P2(cos#) itself has a minimum at 
9 = 90°, the perpendicular coupling at the bandhead is reproduced. These facts lie at the 
origin of the following popular picture that can be found in the literature. States belonging 
to high-j orbitals have doughnut-like shapes with mean angular momenta pointing along 
the symmetry axis of such a density distributions. Therefore, the right-angle orientation 
of the two spin vectors at the bandhead corresponds to a minimal spatial overlap of 
the particle and hole wave-functions. One usually says that this is because the particle- 
hole interaction is repulsive, and tends to minimize the overlap. This reasoning may 
be misleading, because it overlooks the core, whereas the P2(cos#) interaction describes 
precisely the influence of the core polarization, and does not have the sense of a direct 
force. 

With the P2(cos6) interaction at hand, one can address the fundamental question of 
how the energy changes with spin in a shears band. Although the quadrupole correlations 
smooth out the E(I) dependence, the shears mechanism of generating the angular mo- 
mentum is substantially different from that occurring for the collective rotation, and it is 
not obvious whether it should also lead to the rotational-like behavior of E ~ 1(1 + 1). 
In their model study, Macchiavelli and Clark assumed that the P2(cos#) force is the only 
source of the energy change along the band and that the whole angular momentum comes 
from the shears closing, that is, I = j. Under these assumptions, the energy, E, of a level 
belonging to a shears band has the form of P2(cos#), 

„ 3 cos 2 9-1 . . 

E , (2.1) 

whereas from the shears geometry shown in Fig. 12.21 one easily obtains 

Substitution of (|2.2j) to ()2.1|) yields a fourth-order E(I) dependence, and the fourth- 
order term is not only a correction, because the quadratic term has a negative sign. For 
given values of j p and j h , the shears angle, 9, can be extracted from the experimental 
spins according to (|2.2j) and the obtained dependence of E{9) can be compared with the 
form (|2.1j) . Such an analysis was done for 198 > 199 Pb [Mac 98bj and 142 Gd [RzaOlaJ, giving 
a satisfactory agreement with (j2.1j) . The fourth-order dependence translates into non- 
constant moments of inertia. This aspect was exploited in |Mac98c| . where a qualitative 
agreement with the moments of inertia for a band in 197 Pb was found. 

The important role of the quadrupole polarization and the quasi-classical behavior 
of the angular momenta of the high-j nucleons make the magnetic rotation suitable for 
the mean-field description within the cranking model. In this approach, the alignment 
properties of the valence particles and holes are governed by the deformation of the mean 
potential. It turns out that the high-j particles and holes align their mean angular mo- 
mentum vectors, j p and j , on the short and on the long axis of the deformed nucleus, 
respectively, giving the perpendicular coupling at the bandhead. With increasing rota- 
tional frequency, the two blades gradually align towards the axis of rotation, in agreement 
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with the shears picture. Since the valence nucleons provide components of the angular 
momentum on two different axes, the considered system is an example of nuclear rotation 
about a direction tilted with respect to the principal axes of the nucleus. The variant of 
the cranking model that allows for such a rotation is called Tilted- Axis Cranking (TAC), 
and was originally developed by Frauendorf jFra93aj . In fact, the very idea of the shears 
mechanism was first inspired by TAC considerations, and a great majority of experimental 
data on the magnetic bands has been successfully reproduced within this theory. From 
TAC calculations, (3 deformations as small as ~ 0.05 and rarely approaching ~ 0.19 are 
assigned to the shears bands. It seems, however, that smaller deformations are needed 
to reproduce the properties of shears bands than those that come from the energy mini- 
mization. Recently, Frauendorf suggested Fra03aJ that the direct interaction between the 
valence nucleons, absent in the mean field, may be of some importance, or that isovector 
deformations may play some role. Self-consistent methods are necessary to address the 
latter point. The Skyrme-Hartree-Fock (Skyrme-HF) TAC approach developed in this 
work is one of the first theoretical tools appropriate for a study of such effects. 

In general, in generating the angular momentum there is a competition between closing 
of the shears and the collective rotation. Two elements are decisive here: deformation 
and pairing correlations. If deformation is small, there is almost no collective motion and 
the whole inertia comes from the aligning blades. As deformation grows, obviously the 
core rotation becomes more important. But it also turns out that the angular momenta 
of the valence nucleons get more and more rigidly fixed with respect to the nuclear shape, 
retain their perpendicular coupling, and thus do not produce any increase in spin. This 
pushes the balance towards the collective rotation and leads to a change in the total 
moment of inertia. That change is not necessarily drastic, because the J^ 1 moments of 
the shears bands are comparable to those of the well-deformed rotational structures. But 
quantities like B(M1)/B(E2) and J^' / B(E2) are radically influenced, mainly due to 
the increase in B(E2). For a given deformation, inclusion of pairing in the calculations 
will generally reduce the collective inertia. It is argued in Section [Ol that by introducing 
additional mixing of the s.p. states, pairing can also soften the deformation alignment of 
the valence particles and holes and thus make the shears close faster. These arguments 
indicate that the pair correlations favor the shears mechanism over the core rotation as far 
as the generation of the angular momentum is concerned. Macchiavelli and Clark studied 
Mac99aj the shears-rotor competition within a classical model of two blades interacting 



via the P2(cos8) force, coupled to and a spherical rotor. By minimizing the energy at a 
given spin they found the balance between the shears and rotor contributions. By scaling 
the model parameters to the lead region they arrived at a universal estimate that the 
shears mechanism is expected to dominate for e <~ 0.12. For the known magnetic bands, 
it is estimated that the core rotation contributes ~ 15% of the total spin. 

For more information on the magnetic rotation and the shears mechanism, see the 
review articles by Clark and Macchiavelli ClaOOaj and by Frauendorf [FraOlaj . A com- 
pilation by Amita et al. [AmiOlaj lists detailed experimental and theoretical data on all 
the known magnetic dipole bands (May 2001) together with complete bibliography. 
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2.2 Chiral rotation 

The history of chirality began in 1848, when Louis Pasteur noticed that molecules of 
some substances exist in two species, which turn the light polarization plane either left 
or right. Half a century later, lord Kelvin called these forms chiral, from the greek word 
X £L Pi meaning hand, in analogy to the left and right hand. / call any geometrical figure, 
or group of points, chiral, and say it has chirality, if its image in a plane mirror, ideally 
realized, cannot be brought to coincide with itself. - Lord Kelvin, Baltimore Lectures on 
Molecular Dynamics and the Wave Theory of Light, 1904. Thus, a system is chiral if it is 
not symmetric with respect to mirror-reflection, S, in any plane. However, the chirality of 
nuclear rotation does not concern the position vector, but rather the axial pseudo-vector 
of angular momentum. Therefore, it is represented by the operator R T , product of the 
time reversal and rotation through 180°. It is easily verified that R T acts on the angular 
momentum like S acts on the position vector. 

The R T symmetry was introduced into consideration in nuclear structure physics in 
1997, by Frauendorf and Meng [Fra97aJ, who were inspired by a coincidence of two cir- 
cumstances. One was that in 1987 Frisk and Bengtsson demonstrated |Fri87aj that in 
triaxial nuclei mean-field cranking solutions can exist with the angular momentum vec- 
tor, J, having non-zero components, J s , J m , Ji, on the short (s), medium (m), and 
long (I) axes, as illustrated in Fig. 12.31 Frauendorf and Meng pointed out that such an 
orientation of the angular momentum violates the R T symmetry. Indeed, if vectors J s , 
J m , J i form, say, a left-handed set 2 (Fig. 12.3a ). then their images in the mirror reflection, 
R T , will form a right-handed set fFig. 12.3b ). Like the left and right hand, the two systems 
cannot be superimposed by any rotation, and are thus chiral. Notice that the discussed 
most general orientation of J with respect to the triaxial shape also violates the signature 
symmetry, R. 

The second premise of Frauendorf and Meng was the observation, by Petrache et al. 
in 1996 [Pet96al . of an almost degenerate doublet of positive-parity AI = 1 bands in 
i34p r App earance of such pairs can be explained in terms of the breaking of the chiral 
symmetry in the following way. Spectroscopic measurements provide us with eigenvalues 
of observables (spin, energy) that commute with R T . Therefore, quantum states observed 
in the laboratory do not violate R T . At a given spin, the two states of the doublet, |+) 
and |— ), can be interpreted as 

1+) = ^=(\C) + \K)) , |-) = ±(\£) - \K)) , (2.3) 

where \C) and \1Z) are the left-handed and right-handed states represented in Fig. I2.31 
Since states \C) and \TZ) are transformed onto one another by R T , 

R T \C) = \K) , R T \TZ) = \C) , (2.4) 

it is easily verified that |+) and |— ) are, indeed, eigenstates of that operator, 

R T \+) = |+> , #V> = -!->• (2-5) 

2 Of course, which handedness is which is a matter of convention. Here, the right-handed-screw orien- 
tation of space and the s-m-l order of the axes is adopted. 



18 



CHAPTER 2. MAGNETIC AND CHIRAL ROTATIONS 




Figure 2.3: Triaxially deformed nucleus with the angular momentum, J, having non-zero 
components, J s , J m , Ji, on the short (s), medium (m), and long (7) axis. The left panel, 
\£), shows the left-handed orientation of those three vectors, and the right panel, \TZ), 
shows the right-handed orientation. The components J s , J m , J\ originate from the spin 
of the valence particle, j p , from the collective angular momentum, R, and from the spin 
of the valence hole, j h , respectively. 

They are also orthogonal 3 . The nuclear many-body Hamiltonian, H, is invariant under 
R T , and therefore the states \£) and \TZ) have equal mean energies. At the same time, 
the matrix element of H between those states is, in general, non-zero, i.e., 

(£\H\£) = (K\H\K) = E , \(£\H\K)\ = A ^ . (2.6) 

As a consequence, the energies of the laboratory states, |+) and |— ), read 

(+\H\+) = E + A, ( — \H\—) = E — A . (2.7) 

Their energy splitting is, of course, due to the interaction between the left- and right- 
handed states. One can see that doublet bands like those observed in 134 Pr may be 
attributed to the breaking of the chiral symmetry. 

One can thus suppose that the appearance of the doublet bands is connected to the 
existence of triaxially deformed states with the angular momentum vector lying outside 
any principal plane. Frauendorf and Meng pointed out that in odd-odd nuclei around 
i34p r suc j 1 states can arise in the following way. In this region, a configuration is easily 
available, in which one proton particle occupies the lowest substate of the /in/2 orbital, 
and simultaneously one neutron hole is left in the highest /in/2 substate. The former 
drives the nucleus towards elongated shapes, while the latter towards disc-like forms. 
The interplay of these opposite tendencies may yield stable triaxiality. In the triaxially 
deformed potential, the considered particle and hole align their angular momenta, j p and 
j h , along the short and long axis of the nucleus, respectively, providing the J s and J/ 
components of J. As expected from the hydrodynamical irrotational-flow model [Boh75a , 
the moment of inertia with respect to the medium axis is the largest, which energetically 

3 More precisely, they can be always made orthogonal by an appropriate choice of phases of the states 
|£) and \K). 
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favors the collective rotation about this direction. Thus, the component J m comes from 
the collective angular momentum, R. This is illustrated in Fig. 12.31 

Since the first observation in 134 Pr, about 12 similar structures have been found in odd- 
odd nuclei in the A ~ 130 region. They are all attributed to the discussed configuration 
of 7r/i}i/2 v ^Hxit R- ecen tly, first observations of possibly chiral bands in the even-odd 
135 Nd [Zhu03aj and in the even-even 136 Nd |Mer02aj were reported in the same region. 
There, the proposed active configurations are 7r ^ ii/ 2 u ^ii/ 2 an( ^ n ^ii/2 u ^hx/v respectively. 
Peng, Meng, and Zhang predicted theoretically |Pen03a| that another island of chirality, 
associated mainly with the configuration ngl/ 2 vg^j 2 , may exist around A ~ 100. Indeed, 
experimental work is underway for 102_106 Rh iS taOSa, and gives promising results. It seems 
that a pair of bands found in 104 Rh |Koi03b| IVam04a] provides the best known example 
of a chiral doublet, because the chiral partners are situated closest to each other. In this 
case, the irg^h vhn/ 2 configuration is proposed. Latest lifetime measurements in 132 La, by 
Grodner et al. |Gro0 4a . provided the first direct information about the absolute B(E2) 
and B(M1) values in the chiral bands; see Section I""""T1 A complete (January 2004) review 
of experimental and theoretical investigations of the chiral bands in particular nuclides 
is given in Tabs. 12.11 and 12.21 See also Figs. 16.11 and 16.21 for sample level schemes and 
Figs. I6.14H6.151 for spin-energy plots. 

Certainly, an important question about the observed doublet bands is whether the 
chiral rotation is the only possible explanation of their appearance. A more standard 
picture would be that they are in fact four AI = 2 bands corresponding to the two signa- 
ture states of the two odd nucleons, one proton and one neutron. The energy separation 
between the bands would then be just the signature splitting. However, the calculated 
signature splitting turned out to be a few times larger than the observed distance [StaOlb . 
Also the possibility that the yrare partner band is a gamma-vibrational excitation built 
on the yrast band has been ruled out for the same reason. However, recently Pasternak 
|Pas04aj managed to reproduce all the three bands observed in 132 La (see Section 16*70) in 
a very simple classical model that invokes only planar rotation. Yet, some of the model 
parameters were fitted to each band separately. 

Two theoretical tools have been used for the description of the chiral rotation: various 
versions of the Particle- Rotor Model (PRM), originally developed by Bohr and Mottelson 
[BohTSa. and the mean-field Tilted Axis Cranking (TAC) model, which is discussed in 
Section "*"~""1 and used in this work. In the PRM approach, the nucleus is represented by a 
rotating deformed structureless core and a few valence particles interacting with it, and 
possibly among themselves. The valence nucleons are usually treated as pure particles or 
holes, while calculations with pairing appeared only recently |Koi03aj . The core is usually 
taken as a triaxial rigid quantum rotor (Davydov-Fillipov model |Dav58aj ). characterized 
by three moments of inertia, J S1 J~ m , J7j, and described by the Hamiltonian 



H, 





(2.8) 



rot 



^3 s ^^Tm 



where R S} R m , Ri are the intrinsic- frame components of the core angular momentum. In 
all available PRM studies, the moments of inertia, J~k, are calculated from the irrotational- 
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J, 



Ji 



Figure 2.4: Schematic drawing of the trajectory followed by the total angular momentum 
vector, J, along planar and chiral bands. Indices s, m, I refer to the short, medium and 
long intrinsic axis, respectively. Arrows marked as j p and j h denote the initial angular 
momenta of the valence proton particle and neutron hole. Regions of planar, chiral and 
axial rotation are marked. 

flow formula [Boh75aJ, 



where k = 1,2,3 corresponds to the short, medium, and long axis 4 , respectively, and 
B is a mass parameter. The important point is that PRM gives wavefunctions of good 
angular momentum, and thus describes the system in the laboratory frame. Therefore, 
it directly yields the splitting between the bands and the transition probabilities. Due to 
the presence of the rigid core, PRM does not take into account polarization of the nucleus 
by valence particles, nor change of shape induced by rotation. These effects are properly 
included in the mean-field TAC model, especially in its self-consistent version. The mean- 
field wavefunctions, in turn, explicitly break the chiral symmetry, and correspond to the 
\C) and \1Z) states in the intrinsic frame. This makes the physical interpretation more 
straightforward, but neither the splitting nor the transition probabilities can be directly 
calculated. The PRM and TAC methods are to an extent complementary. A more 
detailed comparison of the two models in their assumptions and results can be found in 
|Fra96b[ IFra9 7aj . It is noteworthy that the salient features of the chiral rotation can be 
elucidated within a very simple classical model of two gyroscopes, representing the valence 
particle and hole, coupled to a triaxial rigid body standing for the core. This model is 
discussed in |Olh02bl IUim02al I()lh04a] . and in Section El 

The following general structure of a chiral band could be established from the theo- 
retical studies. Figure E31 shows a schematic trajectory of the mean angular momentum 

4 The irrotational-flow formula, as it is given in Eq. I|2.9[) . and with such an assignment of k to the 
three intrinsic axes, is valid for the 7 range from to 60°, which is used in the present work. 




k = 1,2,3 



(2.9) 
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Figure 2.5: A chiral doublet band for the Kg\j 2 u 9g/2 configuration in the A ~ 100 
region, calculated within the PRM [Pen03aj . a) energy in function of spin, b) spin versus 
rotational frequency calculated as U)(I) = [E(I + 1) — E(I — l)]/2. Solid and dashed 
lines connect the points corresponding to odd and even spins, respectively. The regions 
of planar, chiral and axial rotation are marked. 



vector, J, in the intrinsic frame, and Fig |2.5b gives the I(u) dependence calculated within 
the PRM |Pen03a| for a generic irgl/ 2 v 9gi2 chiral band in the A ~ 100 region. As al- 
ready discussed, the valence particle and hole tend to align their spins, j p and j h , along 
the short and long axes of the triaxial shape, respectively. At low spins, the total J 
comes mostly from those valence nucleons, and remains in the plane spanned by the short 
and long axes (marked as 'planar' in Fig 12.4)1 . This is called planar rotation, in which 
the chiral symmetry is not broken. In this part of the band, the dependence of spin on 
rotational frequency is approximately linear (marked as 'planar' in Fig 12.5b ). with the 
second moment of inertia close to the moments, J s and J\ y for the short and long axes 
(in most studies, 7 ~ 30° is assumed, for which J s m J{). At a certain value of the 
rotational frequency, the spin vector moves out of the s-l plane into one of the octants 
of the intrinsic frame, and the rotation becomes chiral (marked as 'chiral' in Fig 12. 4|) . 
In the current study, an analytical estimate for that critical frequency has been found in 
terms of the classical model; see j()lb02bllOTh04a] and Section El In the chiral regime, 
the spin-frequency plot is also a straight line (marked as 'chiral' in Fig |2.5b ). but is 
now approximately equal to the moment of inertia, J m) associated with the medium axis. 
This results in a kink in Fig. 12.5b around spin 11 %. As the angular momentum vector 
traverses the selected octant further along the band, it aligns more and more with the 
medium axis (marked as 'towards axial' in Fig |2.4|) . Simultaneously, contributions from 
the valence particle and hole become negligibly small as compared to the total spin, and 
one effectively arrives at the axial rotation about the medium axis. In this regime, not 
only the chiral symmetry, but also the signature is unbroken. In view of such a scenario, 
it is stressed in the literature that the chiral rotation is a transient phenomenon, existing 
only in a limited range of angular momentum. 

Figure E31 shows the mean intrinsic-frame trajectory of the angular momentum vector 
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along a chiral band. Classically, each point of that trajectory corresponds to a uniform 
rotation of the nucleus about a fixed direction, without any wobbling. Within the TAC 
model, the trajectory follows the minimum of the energy surface with increasing rotational 
frequency. The two chiral partners can be viewed as two lowest quantum states in the 
potential represented by that surface. In the single planar minimum existing at low 
spins, the two states are interpreted as zero- and one-phonon excitations separated by 
a finite energy interval. They correspond to oscillations of the spin vector about the 
planar equilibrium, which has a classical analogy in non-uniform rotations |Dim02aj . In 
the chiral regime, there is the left and right minimum, of equal depths. If the barrier 
between them is high enough, the two lowest eigenstates are nearly degenerate, and one 
can talk about strong chiral symmetry breaking. For low barriers, the spin vector can 
tunnel between the minima, and the splitting between the band states is non-zero 5 . This 
kind of weak symmetry breaking has been dubbed chiral vibrations [StaOlb . As the 
axial rotation is approached, the weakened breaking of both the chiral symmetry and 
the signature manifests itself, in particular, by the onset of the signature splitting. That 
splitting is absent for the planar and chiral rotation, because the signature is strongly 
broken in those modes. The above features appear clearly in the structure of the chiral 
band, as calculated in the PRM model. This can be seen from Fig. 12 .5[ showing the plots 
of energy versus spin and spin versus rotational frequency. 

The evolution from planar to chiral to axial finds a manifestation in the electromagnetic 
properties of the chiral bands, too. Those were examined in |Fra97a[ IPen03a[ IPen 03b , 
mainly in the frame of the PRM model. For planar rotation, the transitions within 
each partner band are strong, the transitions from the yrare to the yrast band are weak, 
and there are practically no decays in the opposite direction. In the chiral regime, all 
reduced probabilities assume moderate values, and the intraband and interband transition 
strengths are comparable. As the planar rotation is approached, the B(E2) within the 
two partners become equal and slowly increase, while the interband E2 decays disappear. 
The B(M1) show the characteristic odd-even staggering, absent at lower spins. Reference 
|Koi03a| gives a discussion of the staggering in B{Ml)/B(E2) and B(Ml) in /B(Ml) out 
(interband/intraband) in connection with the signature splitting and inversion. 

The appearance of the chiral bands is considered a strong evidence for the existence of 
triaxial deformations; certainly, the chiral geometry cannot occur in an axially symmetric 
nucleus. It is clear from the irrotational-flow formula (J2.9)) that 7 = 30° gives the largest 
moment of inertia with respect to the medium axis {J m = AJ S = AJi) and thus favors the 
aplanar orientation of the angular momentum. The PRM study by Starosta et al. [StaOlat 
ISta02aj showed that the value of 7 practically does not influence the alignments of the 
valence particle and hole on the short and long axis. However, those alignments are only 
well defined for sufficiently large (3 deformation. References [StaOla, Sta02aJ investigated 
the influence of 7 and (3 on the mean value of the so-called orientation operator, a = 
(j p x j h )R, which measures the aplanarity of j p , j , and R, i.e., the degree of chirality, 
in a sense. The average value of a increases with (3 and has a maximum for 7 = 30°. This 
is reflected in the splitting between the chiral partners, which decreases with (3, and is 
minimal for 7 = 30°. However, it is known that the flow of matter in nuclei is not purely 



5 However, in order to actually calculate the chiral splitting in such an approach, one would need the 
mass parameter in addition to the potential surface, and such calculations have not been done. 
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irrotational, i.e., the curl of the velocity field is not exactly zero. Thus, the moments of 
inertia may deviate from the irrotational-flow values (|2.9|) assumed in the PRM studies, 
especially if pairing is weak, and the dependence of the discussed observables on 7 may 
not be so straightforward. 

As far as the mean-field methods are concerned, certainly taking into account the 
interaction between the left- and right-handed minima would be desirable, preferably 
in conjunction with the projection onto good chirality before variation. This would be 
the first tool capable of calculating the chiral splitting and the transition probabilities 
in a fully microscopic way. In the PRM domain, efforts are being undertaken to include 
pairing in the calculations |Koi03aj , to use 7- soft cores, and to take into account the direct 
residual interaction between the valence nucleons |Rai03a| IRai03bj . Such an interaction 
is supposed to be responsible for the stabilization of the chiral geometry in nuclei at the 
limits of the islands of chirality. Experimentally, lifetime measurements would certainly 
be of great use. 
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Table 2.1: Experimental and theoretical investigations of candidate chiral rotational 
bands; this Table is continued in Table 12.21 Acronyms in column Theory stand for calcu- 
lations in the following models. TAC: phenomenological Tilted-Axis Cranking model of 
Frauendorf et al; PRM: Particle- Rotor Model used by Peng et al. [Pen03a, Pen03E] or 
its different version used by Hartley et al. [HarOlaj : CPH: Core- Particle-Hole model, a 
version of PRM developed by Starosta et al] CQP: Core-Quasi-Particle model, a version 
of CPH that includes pairing; SHF: Skyrme-Hartree-Fock results of the present work; see 
Tab 16 .21 for details. The symbol v in the column Pairing means that pairing was included 
only for neutrons. 
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Table 2.2: Continuation of Table l2~Tl 
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CHAPTER 2. MAGNETIC AND CHIRAL ROTATIONS 



Chapter 3 
Theoretical tools 



In this Chapter, some general information about the spontaneous symmetry breaking in 
the mean field is recalled, mainly in the context of rotational bands. While the Tilted- 
Axis Cranking (TAC) method has already been widely discussed and used in its version 
based on phenomenological potentials, this work introduces one of its first applications 
to the self-consistent Hartree-Fock (HF) calculations. Therefore, the TAC basics are 
recapitulated with particular emphasis on the differences between the phenomenological 
and self-consistent implementations. It is also described how the values of the time-odd 
coupling constants of the HF energy functional are chosen in the present calculations 
to investigate the role of the time-odd nucleonic densities and fields in the shears and 
chiral solutions. After a concise review of other computer codes capable of performing 
symmetry-unrestricted mean-field calculations, the program hfodd, used in this work, is 
presented. 

3.1 Spontaneous symmetry breaking 

This Section summarizes the basics of the HF method, and explains how the spontaneous 
symmetry breaking intervenes in the mean-field approach. 

The variational HF method consists, in its standard version, in minimizing the expec- 
tation value of the many-body Hamiltonian, 



in the trial class of Slater determinants. Here, T is the kinetic energy, and V is the 
effective two-body interaction. A Slater determinant, |\&), is uniquely characterized by 
its hermitian, projective density matrix, 



H =T+V , 



(3.1) 




(3.2) 
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By evaluating the derivative one obtains 

h[p]=f + t[p} , (3.4) 

where T[p] is the mean potential generated by averaging the two-body interaction with 
the density distribution, 

T[PW = ^ VftPvaPaP ■ (3-5) 

a/3 

The HF equations state that a local energy minimum is attained for such a density matrix, 
which commutes with the s.p. Hamiltonian it induces, 

p,h\p]]=0. (3.6) 

The last equation tells us, in particular, that the density matrix and the s.p. Hamiltonian 
have common eigenstates. Therefore, the Slater determinant corresponding to the HF 
solution represents a set of particles moving in the mean potential T[p\. Although it 
is known that a Slater determinant provides a poor approximation to the nuclear stat, 
the obtained approximation to the s.p. density, p, is much better. Therefore, like in 
the Density Functional Theory, of which the HF method is a particular reliable 
physical result is contained in the density. 

The nuclear many-body Hamiltonian, H, is invariant under a number of symmetry 
operations, S, 

[H, S] = . (3.7) 

Normally, these conserved symmetries include the translation, rotation, plane reflection, 
time reversal, and particle-number symmetries. This list can be supplemented with the 
approximate isospin symmetry, which is broken mostly by the Coulomb interaction, and 
by the Galilean invariance, which is only a symmetry of the interaction V. In different 
physical conditions, e.g., when the system is subjected to an external electromagnetic field, 
some of those symmetries are unconserved. The nuclear states belong to irreducible rep- 
resentations of the group formed by the conserved symmetries of the Hamiltonian H. In 
general, those representations may be multi-dimensional, like for the angular-momentum 
eigenstates with L > 0. Such states are said to be covariant with the group operations. 
States that are bases of one-dimensional representations, like the angular-momentum 
eigenstates with L — 0, are called invariant 1 . In the present work, only conservation 
of single dichotomic symmetries, like plane reflections, is analyzed. A dichotomic sym- 
metry is a two-element group containing the given symmetry operation and the neutral 
element. Such a group has only two irreducible representations, both one-dimensional, 
the symmetric and the anti-symmetric one. 

In the variational HF approach, one does not calculate exact eigenstates of H, but 
looks for their approximation among possibly simple states, the Slater determinants. If, 
for instance, there are strong quadrupole correlations in the exact solution, the only way 
the simple determinant can account for them may be by a static quadrupole deformation, 
which violates the rotational invariance. Thus, in general, the mean-field solution does 
not possess the symmetries of the many-body Hamiltonian, which means that the density 



1 Also for one-dimensional representations different than the totally symmetric one. 
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matrix, p, corresponding to the minimum energy, does not commute with the symmetry 
operator, S, 

[p, S}^0; (3.8) 

then one says that the symmetry S is spontaneously broken in the HF solution. But it 
may also remain unbroken, which depends on the constituents of the system and on the 
interaction V. It is clear that symmetries broken in the mean field are reminiscent of the 
presence of certain correlations in the exact solution. Although the quantum mechanics 
formally requires a proper symmetry behavior, one may note that the mean-field descrip- 
tion eventually becomes exact in the limit of large systems, like macroscopic crystals or 
even molecules. Thus, it is also sufficient for the description of certain properties of nuclei. 
If information of a more quantal character is required, like the excitation spectrum, one 
may need to restore the broken symmetries, e.g., by projecting the mean-field solutions 
onto good angular momentum. See |Rin00a| for a more thorough discussion. 

An important result about the symmetry breaking in the HF approach is provided by 
the theorem of self-consistent symmetries [RinOOaj , which requires that if S is a symmetry 
of the many-body Hamiltonian, H, then 

S+h[p]S = h[S + P S] . (3.9) 

In particular, if p is invariant under S then so is h. On the other hand, if h commutes 
with S then its s.p. eigenstates can be chosen as eigenstates of S, which means that also 
p is invariant under S. Consequently, the s.p. density and the s.p. Hamiltonian have the 
same symmetries in a self-consistent solution. In the present study, the condition of the 
theorem are always satisfied, because the Skyrme interaction is rotationally invariant, and 
only subgroups of the rotation group are considered here. 

The exact solutions are often referred to as laboratory states, which are, in general, 
covariant with the symmetry group of the many-body Hamiltonian. One says that spon- 
taneous symmetry breaking occurs in the intrinsic frame, which actually alludes to the 
mean-field solutions. In Section l2~2~l the states |+) and |— ) of a chiral doublet are the lab- 
oratory states, and the left- and right-handed states \C) and \1Z) are mean- field solutions, 
while the passage from the latter to the former is the projection onto good chirality. 

3.2 Symmetries and rotational bands 

Relation between the symmetries of the nucleus and the structure of rotational bands is 
presented in this Section. The symmetry group Dj h , which plays a crutial role in the 
present analysis, is defined. 

For all quantum systems whose rotational excitations are observed, like molecules or 
nuclei, one can find universal connections between the symmetries in the intrinsic frame 
and the sequence of spin and parity of levels within a rotational band. This can be done 
by checking what values of angular momentum and parity can be projected out of a mean- 
field solution with given symmetries unbroken. Frauendorf |Fra01a| IFraOlbj summarized 
the results of such an analysis in a table which is reproduced in Table l3~Tl It can be seen 
that three symmetry operations are relevant for the classification of the bands, namely, 
the space inversion or parity, P (reflection in a point), the signature, R (rotation through 
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Table 3.1: Symmetries of the nucleus in the intrinsic frame and sequences of spins and 
parities in the corresponding rotational bands. The parity (P), signature (R z ), and chiral 
symmetry (Ry) operators are defined in the text. Symbols U and B stand for unbroken 
and broken, respectively; an operator as entry means that this operator gives the same 
result as the one given in the header of the column. 7+ means that there are two degenerate 
states of opposite parity. The 2 means that there are two degenerate states with the same 
7+. In rows I-V, the parity of the band can be positive or negative, although the latter 
possibility is not explicitly indicated. From jFraOlal Fr aOlbj . 
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R y 


Level sequence 
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7+,(7 + 2)+, 


(/ + 4)+ 


TT 
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T T 

u 
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XV 


B 


B 


B 


2/+,2(/ + l) ± 
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180°), and the chiral symmetry, R T (product of the signature and the time reversal). 
Not all of the 15 mathematically conceivable kinds of bands listed in the Table have 
been experimentally found, and the interesting question whether they all appear in nuclei 
remains open. 

Together with the important time-reversal operation, the symmetries appearing in 
Table 13.11 generate a group called D\ h [DobOOb , which is of crucial importance in the 
present analysis. The group comprises 2 : 

• Three signature operations, Ri = exp(—inJi), i = x,y,z, which are rotations 
through 180° about the three cartesian axes. Here, Jj denotes components of the 
angular momentum operator. 

• Space invertion, P, sometimes called parity, changing r into —r. 

• Three simplex operations, Si = Pexp(— iirJi), which are reflections in the planes 
perpendicular to the i = x,y,z axes. 

2 The identity element, 1, the change of sign, —1, and the products of —1 with other elements of the 
group are not mentioned. 
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• Antilinear time-reversal operator, T = exp(— ms y )K, where s y is the y component 
of the total intrinsic spin operator, and K is the complex conjugation in the space 
representation. 

• T-signatures, Rj = TRi, T-parity, P T = TP, and T-simplexes, Sf = TSi, where 
i = x } y } z. 

The T-signatures are the chiral operators discussed in Section 

When enumerating the Dj h symmetries it is necessary to introduce names of the three 
cartesian axes. But obviously, the question in the mean-field calculations is not, e.g., 
whether the signature with respect to some particular a priori given axis is broken or not. 
One says that the signature is broken uniquely if there is no axis % such that Ri is the 
symmetry of the density. Otherwise, the signature is not broken. The same pertains to the 
simplex, T-signature, and T-simplex symmetries. In the following, the terms broken and 
unbroken are always understood in this sense, when referring to the signature, simplex, 
T-signature, and T-simplex symmetries, unless a particular axis is specified. 



To obtain HF solutions describing states with higher angular momentum, particularly the 
levels of a rotational band, one usually imposes a linear constraint on the total angular 
momentum, and minimizes the mean value of the quantity called Routhian, 



where the Lagrange multiplier uj is called rotational frequency, and J is the angular- 
momentum operator. This approach is called cranking approximation. Since H is rota- 
tionally invariant, the solutions obtained for the same length, but different directions of 
uj differ only by their orientation in space. Therefore, only the length, u, has a physical 
meaning. The Kerman-Onishi theorem |Ker81aj states that, at the minimum, the total 
angular momentum vector, J = (J), is parallel to u>. Indeed, this configuration obviously 
minimizes the mean value of the cranking term, — ojJ, while H is rotationally invariant. 
It comes from the variation of the many-body Routhian, H', that in the cranking model 
the s.p. Hamiltonian, h, is replaced by the s.p. Routhian, 



both in the HF equations ()3.6|) and in the theorem of self-consistent symmetries (j3.9|) . 

If a given symmetry is not broken in a mean-field solution, expectation values of observ- 
ables represented by operators odd under this symmetry vanish. Here, we are interested 
in properties of the D\ h symmetries with respect to the mean angular momentum vector 
and electric quadrupole moment of the mass distribution, because they will play a fun- 
damental role in the search for broken signature and chiral symmetries. These properties 
can be summarized as (see |Dob00b| for a detailed derivation): 

• If the time reversal, T, is not broken, all components of the angular momentum 
vanish, while the quadrupole moments are not affected. 



3.3 Tilted- Axis Cranking 



H' = H — u)J , 



(3.10) 



h! = h — uj J , 
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• Unbroken parity, P, does not constrain neither the angular momentum nor the 
quadrupole moments. 

• If Ri or Si is unbroken, the angular momentum can have a non-zero component only 
on the i-th axis, which in turn is a principal axis of the quadrupole tensor. 

• If Rf or Sf is unbroken, the axis i is a principal axis of the quadrupole tensor, while 
the angular-momentum vector is confined to the plane perpendicular to that axis. 

These conclusions are valid, of course, for any arbitrarily oriented axis %, not only the x, 
y, z axes of some particular frame. It follows from the third bullet that if the angular 
momentum has non-zero components on at least two principal axes of the quadrupole 
tensor, then there does not exist any axis i such that Ri or Si be an unbroken symmetry. 
The last point tells us that if the angular momentum has non-zero components on all the 
three principal axes of the quadrupole tensor, then there is no such axis i that Rf or Sf 
be unbroken. 

It is clear from the above properties of the Df h symmetries that whenever there exist 
an axis % such that either R iy Si, Rf, or Sf is an unbroken symmetry, then i is a principal 
axis of the quadrupole tensor. So, only the principal axes and planes of this tensor 
are candidates to being symmetry axes and symmetry planes. It is, therefore, useful to 
consider the intrinsic frame of reference, defined by the principal axes of the quadrupole 
moment. To denote the axes of the intrinsic frame, the symbols s, m, I, are used, in view 
of the short, medium and long axes of a triaxially deformed shape. Let us also introduce 
the program frame, which is defined by the axes used for solving the mean-field equations, 
e.g., in a computer code 3 . They are denoted as x, y, z. In general, the axes of the two 
frames do not coincide. 

If one wants to check by using the mean-field methods whether a given symmetry 
is spontaneously broken or not in the considered system, one looks for the minimum of 
the Routhian among trial solutions both breaking and not breaking the symmetry in 
question. Then one checks, to which of those two categories the minimum belongs. These 
are the symmetry-unrestricted calculations. One can also impose a given symmetry, which 
amounts to searching for the Routhian minimum only among solutions that do not break 
that symmetry. When performing calculations with imposed Df h symmetries one can, 
without loss of generality, impose the symmetries associated with the axes of the program 
frame, e.g., R x , which is technically much simpler. In fact, imposing either Ri, Si, Rf or 
Sf, where % — x,y, z, causes that one of the intrinsic axes coincides with the program i-th 
axis. In the following, the term imposed always has this meaning. 

If some symmetries are imposed in the calculations, notably the cranking term, —a; J, 
must conform to them. This means that all the above restrictions concerning the mean 
angular momentum, J, equally touch the cranking frequency, u). For example, solutions 
invariant under the time reversal may only result if no cranking is applied. These remarks 
allow us to introduce the following widely used classification of calculations with imposed 
symmetries. Whenever a signature Ri or simplex Si symmetry is imposed, one can talk 

3 It would be unfortunate to call it laboratory frame, because in the common interpretation of the 
cranking model one would say that the program frame rotates in the laboratory system. The program 
and intrinsic frames do not move with respect to each other, they differ only by the orientation of their 
axes. 
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about one-dimensional cranking, because the applied rotational frequency, u>, together 
with all resulting (single-particle (s.p.) and total) angular momenta vectors, point along 
the unique, i-th axis. A name Principal-Axis Cranking (PAC) is also adopted, because 
the i-th axis is then a principal axis of the mass distribution (its quadrupole moment). 
In PAC calculations, the Kerman-Onishi condition (J \\ uj) is satisfied identically. In 
the Tilted- Axis Cranking (TAC), all signature and simplex symmetries are broken, and 
therefore the frequency, u), and the mean angular momentum, J, vectors can become 
tilted with respect to the principal axes. If either a T-signature or T-simplex is still 
imposed, these vectors remain in a single plain, and hence the name of planar or two- 
dimensional cranking (2D) is used. If not, uj and the resulting angular momentum, J, 
can take any direction, and one talks about three-dimensional cranking (3D). In both the 
2D and 3D cases, still the parity can be either broken or not. 

To solve the HF equations, one often uses the iterative method [RinOOa . It follows 
from the theorem of self-consistent symmetries (see Section 13. 1J) that if a symmetry is not 
broken in a certain iteration, it will remain so in all subsequent ones, and consequently 
in the final result. Therefore, to impose a symmetry it is enough to start the iterations 
from an h even under this symmetry and also to apply the cranking frequency w in a 
proper direction, as discussed above 4 . To render the calculations symmetry-unrestricted, 
either the initial h or the cranking term must violate the considered symmetry. If the 
symmetry is not dynamically broken in the physical solution, it will recover in the course 
of the iterations. In the present calculations, a Z)J h -symmetric Nilsson |Nil55aj potential 
is always used as the first approximation to h, unless the iterations are restarted from 
some previously obtained solution. Therefore, the Dj h symmetries can be broken only 
by the cranking term. If one wants to perform symmetry-unrestricted calculations for 
a non-rotating state, one may nevertheless apply a small cranking frequency in the first 
iteration, only to break the D\ h symmetries, and then continue with no cranking. 

To relate the discussion to the magnetic and chiral rotation, after having performed 
a symmetry-unrestricted calculation one should verify whether the signature or chiral 
symmetry is broken by considering the mean angular momentum vector in the intrinsic 
frame. If it has non-zero components on more than one intrinsic axis, then the signature 
symmetry is broken. If all three components are non-zero, also the chiral symmetry 
is broken. However, to make one sure that these symmetries are unbroken is not so 
straightforward, because there may a priori exist other factors violating them, different 
that the angular momentum coupling, e.g., exotic deformations. A detection of unbroken 
symmetries based on the values of the variational parameters (e.g., elements of the density 
matrix) would be doable, but difficult. Another way is to perform a separate calculation 
imposing the symmetry in question and then to check, whether the results are equivalent, 
i.e., whether they differ at most by orientation in space. This can be verified by looking at 
the total energy, contributions thereto from the different terms of the Skyrme functional, 
density integrals, various observables, and s.p. energies. 

A phenomenological approximation to the HF approach consists in replacing f with 
a model potential, like the Nilsson |Nil55aj or Woods-Saxon |Woo54a] one, which de- 
pends on several parameters, mainly multipole deformations, a\^, of the nuclear surface. 



4 It is still possible that an initially unbroken symmetry breaks due to numerical noise. To avoid such a 
noise, in practical calculations the terms of h' odd under the considered symmetry are simply disregarded. 
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Then, the expectation value of the Routhian h! (corrected for the liquid drop energy via 
the Strutinsky renormalization [RinOOa ) is minimized with respect to parameters of the 
potential. 

It is well understood that only the relative orientation of the angular momentum vector 
and the intrinsic axes carries a physical information. Their position with respect to the 
program frame is meaningless and can be arbitrary. One can take advantage thereof 
in the phenomenological method, where the orientation of the mean potential is under 
control via the variational parameters, like the multipole deformations «a m of the nuclear 
surface, and can be fixed so that the tensor be diagonal in the program frame simply 
by setting an = Im CK22 = 0. With this condition satisfied, the intrinsic frame coincides 
with the program frame, facilitating many computational details. The minimization of 
the expectation value of the Routhian h! at a given magnitude of u is then performed by 
varying the direction of u) and all ax^ except «2i and Im «22- The position of uj with 
respect to the coinciding intrinsic and program frames is often described in the spherical 
coordinates, if, called tilt angles. It is in fact this procedure that is usually understood 
under the term TAC and that has been described many times in the literature, see e.g. 
the seminal papers by Frauendorf et al. [Fra91a, Fra92a, Fra93a, Fra 97c| IFraOOaj and 
|Dim00hj . 

In the HF method, however, the Euler angles defining the orientation of the intrinsic 
axes are not free variational parameters, but complicated functions of the densities, that 
change from one iteration to another according to the HF recipe. Thus, the orientation 
cannot be easily controlled. The only feasible way to keep the nucleus in the axes of the 
program frame is by imposing dynamical constraints on the off-diagonal components of 
the quadrupole tensor and requiring that they vanish. Then, one can vary the cranking 
frequency vector like in the phenomenological approach. It is actually the only way to 
proceed if the energy dependence on the tilt angles, <p, is sought. It has been checked, 
however, that multipole constraints strong enough to confine the nucleus may easily lead 
to divergenciesm, at least in the present calculations. If only the minimum is of interest, 
a more natural and faster way is to fix the rotational frequency vector with respect to the 
program frame and let the mean potential reorient and conform to it self-consistently in 
the course of the iterations. It is but this procedure that corresponds to the minimization 
of the given Routhian (j3.10|) . The direction of the rotational frequency vector in the 
program frame no more plays any role, while the solution now does become tilted with 
respect to that frame. This method has been adopted in the present work. 

Some quantities, like mean angular momenta and multipole moments, carry a clear 
information only when expressed in the intrinsic frame of the nucleus. But they are first 
calculated in the program frame. Since the two no more coincide, it is necessary to find the 
principal axes of the nuclear shape by diagonalization of the electric quadrupole tensor 
and to transform the considered quantities into the corresponding frame by use of the 
Wigner matrices. 

To examine deformed nuclei one usually choses deformed bases. However, deformation 
of the basis improves the description only if it suites the shape of the nucleus; otherwise 
it may spoil the result. For instance, deformation of the Harmonic-Oscillator (HO) basis 
amounts to stretching the cartesian coordinates associated with the axes of the program 
frame, and to taking different numbers of quanta in each direction. This is not appropriate 
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for representing the TAC solutions, whose principal axes are arbitrarily tilted with respect 
to the program frame. A deformed basis would drive such a solution like a force to 
a certain orientation different from the self-consistent one, determined by the direction 
of the cranking vector. In particular, the Kerman-Onishi condition would be violated. 
Therefore, in self-consistent TAC calculations, the s.p. space spanned by the basis should 
be rotationally invariant. It means that the basis should have the three characteristic 
frequencies equal and comprise only entire HO shells. To obtain a reasonable description 
for deformed nuclei, the only way is to use large bases. 

It follows from the theorem of self-consistent symmetries that the s.p. Routhian and 
the density matrix are invariant under the same symmetry operations; see Section 13.11 
However, one should be aware of one exception to that rule, that concerns the time- reversal 
symmetry in case the time-odd terms of the mean field are artificially set to zero in the 
calculations. This simplification is always the case in the phenomenological mean-field 
approach, where time-odd terms simply do not appear in the model potential. With the 
time-odd fields cut off (and with no cranking), the Routhian is explicitly time-even, while 
the density will still acquire time-odd components whenever there is an unpaired particle; 
see Appendix E] The spectrum of the time-even Routhian exhibits the two-fold Kramers 
degeneracy, and the state of the unpaired particle can be any unitary combination of the 
states in the Kramers pair; the physicist must take the decision. The main point of this 
discussion is that this freedom of choice influences the time-odd densities, while the time- 
even ones remain undisturbed. This lemma is proven in Appendix El Together with the 
time-odd densities, also the time-odd observables, notably the mean angular momentum, 
depend on that choice. The time-even quantities, like the quadrupole moments, are 
uniquely determined, because they depend only on time-even densities. These conclusions 
are valid for both the phenomenological and the self-consistent implementation of the 
mean-field method. However, in the self-consistent approach, in each iteration there is 
a feedback from the density to the mean field. One may doubt, therefore, whether the 
choice of state for the unpaired particle alters the generated mean field. This is not the 
case, because the time-even fields originate only from the time-even density, which are 
not disturbed, while the time-odd fields are suppressed. 



One important point in the investigations of the magnetic and chiral rotation is to examine 
the role of various time-odd observables, like current and magnetic-moment distributions. 
Their impact on the results is governed by the time-odd coupling constants of the HF 
energy functional. After introducing necessary definitions of the local densities, energy 
functional, mean fields, and local gauge invariance, this Section explains how the values 
of the time-odd coupling constants are set in the present calculations. 

In the HF approximation, the total energy is, in general, a functional of the non-local 
one-body density, 



where is the many-body wavefunction, r and a are the position and spin variables, 
and a = n, p for neutrons and protons, respectively. The density p a (rcr,r'cr') can be 
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written as a sum of the scalar, p a (r, r'), and vector, s a (r, r'), terms, 

p a (r, r') = ^ Pa(ra, rV) , s a (r, r') = ^ p a (ra, r'a')<T a , a , (3.13) 

U (T(j' 

as 

p a (ra, r'a') = \p a (r, r')5 aa > + |s a (r, r') ■ a aa i . (3.14) 

The HF equations are usually solved within the Local-Density Approximation (LDA), in 
which the total energy depends only on local (r = r') densities, and on their derivatives at 
r = r' up to the second order. One has to consider the following real nucleonic densities 
|Eng75a| : 



(3.15) 



Particle: 


PaS 




= Pa(r,r) , 


Spin: 






= s a (r,r) , 


Kinetic: 






= [V- V'p a (r,r')] r=r , , 


Vector kinetic: 


T a \ 




= [V-V's a (r,r')] r=r , , 


Momentum: 






= ^[(V-V')p Q (r,r')] r=r , , 


Spin-current: 




» 





of which p a , r a , and J a = J^u,a are time-even, and s a , T a , and j a are time-odd. For 
each of the local densities, one defines the isoscalar and isovector density, e.g., 

P0 = Pn + Pp Pl = Pn-Pp (3.16) 

The isoscalar and isovector terms are denoted by subscript t — 0, 1, respectively. 

The interaction part of the HF LDA energy functional can be expressed as a space 
integral, 

£ = J2 [ d 3 r (T-tr n + n° t dd ) , (3.17) 
t J 

of energy densities, 7ij vcn and 7i° dd , which depend on the time-even and time-odd nucle- 
onic densities, respectively. They are defined as 



nr n (r) = c p tP i + c t Ap Pt Apt + cj Pt r t + cCj\ + C? J Pt V ■ J t , 

Hf d (r) = C°s* + C t As s t -As t + C?s t -T t + Cti 2 t +c73s t -(Vxj 



(3.18) 



t) i 



where the square of the tensor density is defined as J \ = ^2^ u J^ U j, and its vector part, 
J t , is defined as J\ t = Ylav e ^nv^^,t- The coupling constants C t appearing in TCf en and 
7i° Ad are conventionally called time-even and time-odd, respectively. They are functions 
of the parameters of the Skyrme force [Dob95aj, and may depend on density if the force is 
density-dependent. Note, however, that in the spirit of the Density- Functional Theory the 
energy functional does not need to come from any interaction, and the coupling constants 
Ct play the role of fundamental parameters. 

Upon variation of the interaction energy, £, with respect to local densities, one obtains 
the neutron and proton Skyrme-HF s.p. mean potentials (see Section l3~T|) . 

■p -peven , -podd , -peven , piodd p peven , podd peven podd /o m\ 

J-n — J-o +i +i l +i l L P~ L +i _i l _i l J [6.L\)) 
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where the time-even and time-odd fields, T f j ven and T° dd , read 

peven = _\7 • [M t V] + U t + ^ • t + t ' VeP 

r odd = _ v . [(o .. Ct)v] + (T . St + j. (v . Jt + J / v) 

They depend on the nucleonic densities through the potential functions, 



(3.20) 



U t (r 

M t {r 
Ct{r 
*Bt{r 
It(r 



2C p tPt + 2Cf p Ap t + ClT t + C t VJ V -J t + U t 
2Cfs t + 2Ct s ^s t + C?T t + Cf- 7 V x j t , 

C t s t , 

2Cfj t - C t VJ V Pt , 
2C|j t + C^V x ^ . 



(3.21) 



The tensor gradient operators in (|3.2U|) and ()3.21|) are defined as (Vcr)^ = V ^Oy and 
V Mi , = ^x^^xV \. The term U' t represents the rearrangement terms resulting from the 
density dependence of the coupling constants, Ct- In standard parametrizations of the 
Skyrme interaction, which are used in the present work, only Cf and Cf depend on the 
isoscalar particle density, po, and then 

V {r) = Y,(^A + %?4) U[(r) = 0. (3.22) 



t' 



d P o rt dp, 



If the energy functional is derived from the Skyrme interaction, there are unique 
relations between the time-even and time-odd coupling constants, Ct- In particular, 

c? = -qr, c J t =-cJ, c7 3 = +c7 J . (3.23) 

However, it turns out [Dob95aJ that the three relations ()3.23|) have a deeper sense, namely, 
they are necessary and sufficient conditions for the invariance of the energy functional 
under a local gauge transformation of the A-bodj wave-function, |\&), 



*0 = exp iV^) I*) , (3.24) 




where 4>(r) is an arbitrary real function of position. In particular, for <p = p ■ r /h, the 
gauge transformation ()3.24j) corresponds to a Galilean boost with velocity v = p/m, 
where m is the particle mass. 

In the present calculations, the SLy4 |Cha97aj and SkM* [Bar82aJ parametrizations of 
the Skyrme force were used. Like in the original fits, the time-even coupling constant C( 
and the time-odd constant Cj were always set to zero. This does not hamper the local 
gauge invariance. Apart from C/, all other time-even coupling constants were taken as 
they come from the parameters of the Skyrme force. To examine the role of the time-odd 
densities and fields, attempts were undertaken to obtain each solution with three sets of 
values of the time-odd coupling constants. In each case, values of selected constants were 
taken as they come from the Skyrme force, while those of the remaining ones were set to 
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Table 3.2: Three sets of values of the time-odd coupling constants used in the present 
calculations. For each set, values of the listed constants were taken as they come from 
the Skyrme force, while the remaining constants were set to zero. 



Set symbol 


Coupling 


constants 


N 






G 




d cV 


T 







zero. The constants selected in each set are listed in Table 13.21 The density-dependent 
and independent components of C* were suppressed or not simultaneously. These three 
types of calculations are denoted throughout the text as N, G, and T. In the first 
case, N, the mean field contained only time-even components, like the phenomenological 
mean potentials, and the local gauge invariance was violated. In the set G, only those 
time-odd coupling constants were not set to zero, that must be related to the time-even 
coupling constants according to (j3.23|) for the energy functional to be invariant under 
the local gauge transformations (|3.24|) . Finally, in the set T all the time-odd coupling 
constants were taken as they come from the Skyrme force, with the exception of Cj. 
In the following Chapters, setting or not selected time-odd coupling constants to zero is 
referred to as excluding or including the corresponding time-odd terms of the mean field, 
as given by Eqs. (l3~27H l3~2Tl) . 

The results obtained for each of the sets, N, G, T, are compared for particular bands 
in Chapters El and El A general conclusion is that the inclusion of the time-odd fields does 
not change the solutions qualitatively. Switching them on and off influences mostly the 
moments of inertia. Additionally, convergence problems arise at low cranking frequencies 
for the G, and particularly the T set, and no solutions could be obtained in several cases. 



3.5 Symmetry-unrestricted mean-field codes 

This Section gives a review of computer codes that can perform symmetry-unrestricted 
mean-field calculations, particularly within the TAC model, and that existed prior to the 
present study. 

Imposed symmetries in mean-field calculations allow to reduce the sizes of matrices 
or of the space lattice, and thus to gain on the computer memory and CPU time. Before 
the interest in magnetic bands emerged at the beginning of the 1990's, the signature or 
simplex symmetries were inherently imposed in the majority of mean-field codes. One 
of the first numerical TAC calculations were carried out in 1987 by Frisk and Bengts- 
son |Fri87aj by using the Nilsson potential Nil55a with pairing. The phenomenologi- 



cal TAC code originally written by Frauendorf Fra93aJ IDimOObj is now widely used for 



description of experimental data on the magnetic and chiral bands. It allows for in- 
cluding monopole pairing correlations with quasiparticle excitations. Recently, Dudek 
and Schunck I )ud() la! have written a relativistic Woods-Saxon TAC program, and first 
results are about to appear soon. One of the first symmetry-unrestricted calculations 
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within realistic self-consistent methods were performed in 1991 by Umar et al. |Uma 91a . 
Without cranking and pairing, they solved the Skyrme-HF equations in the space rep- 
resentation, by using the Spline Collocation Method. About the year 2000, Yamagami 
and Matsuyanagi developed a symmetry-unrestricted Skyrme-HF program working on a 
Cartesian mesh with the imaginary-time method VaiiiOOal lYamOla] . They include ei- 
ther only cranking or only pairing with the density-dependent delta force and without 
quasiparticle excitations. Calculations of Madokoro et al. for a magnetic band in 84 Rb 
MadOOaJ were done in 2000 by use of a Relativistic Mean Field TAC code imposing only 
the parity, working in the harmonic-oscillator basis and not including pairing. For the 
purpose of the present work, the author developed a Skyrme-HF TAC program, hfodd 
version (v2.05c), that works in the harmonic-oscillator basis and allows for breaking of all 
the Dj h symmetries. Currently, the pair correlations cannot be taken into account in the 
TAC mode. This program is described in the next Section. 

3.6 The program HFODD (v2.05c) 

In order to perform self-consistent calculations for the magnetic and chiral rotation, the 
author of this work developed a Skyrme-HF computer code, hfodd version (v2.05c). In 
the previous versions, (vl.60r) and (vl.75r) [DobOOaJ, the simplex symmetry was imposed, 
thus not allowing for TAC calculations. Removing this restriction was the principal im- 
provement introduced during the present work, and led to the version (v2.05c), which 
was used for obtaining results presented here. A brushed-up version, (v2.08i), was pub- 
lished in [Dob04a . This Section gives an overview of the program's main features, with 
particular emphasis on those that have been introduced for the purpose of the present 
investigations. 

The code hfodd solves the nuclear HF equations for the Skyrme effective interaction 
by using the iterative method. The s.p. wavefunctions are expanded onto the deformed 
Cartesian Harmonic-Oscillator (HO) basis. From the s.p. wavefunctions, the HF local 
densities (j3.15|) . and then the mean potentials ()3.20|) are calculated on a Cartesian mesh. 
By using the Gauss- Hermite quadrature, matrix elements of those potentials between the 
HO states are calculated. The matrix of the s.p. Hamiltonian, obtained in this way, is 
then diagonalized, to find the s.p. wavefunctions in the next iteration. The Nilsson po- 
tential |Nil55aj is implemented as a starting point. All densities f!3.15|) appearing in the 
Skyrme functional are calculated, and all time-even and time-odd fields ()3.20|) are taken 
into account. Several parametrizations of the force are available and can be modified. 
The generalized spin-orbit terms Rei95aJ can be included in the Skyrme functional. Ap- 
propriate time-odd coupling constants of the energy functional can be determined either 
from the Skyrme force or from the Landau parameters |Ben02aj . For the center-of-mass 
motion, either the standard one-body correction before variation or a two-body correc- 
tion after variation can be used. Apart from such basic quantities like s.p. and total 
energies, the code calculates the angular momentum vectors, electric, magnetic, and sur- 
face multipole moments, mean-square radii and Bohr deformation parameters. Linear and 
quadratic constraints on spin as well as quadratic constraints on the electric and surface 
moments are implemented. With the simplex imposed as an unbroken symmetry, the 
program can also work in the HFB mode, with either the fixed-gap or density-dependent 
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delta pairing. The Fermi level is found automatically by using the equivalent-spectrum 
method [Dob84aJ. The quasi-particle blocking is not yet programmed. 

The user can choose the pattern of imposed symmetries by independently imposing 
or not the S y , R y , P, T, S'J, Sy, and Sj symmetries. Keeping in mind that the unbroken 
symmetries must form a group, 34 combinations listed in Table l3~3~l are permitted. Some 
of them differ only by an exchange of the Cartesian axes associated with the signature 
and simplex symmetries, and are thus physically equivalent. Altogether, 19 out of the 28 
non-equivalent subgroups 5 of Dj h can be selected as imposed symmetries. Imposing the 
missing 9 subgroups can easily be implemented provided there is a physics motivation to 
study such new cases. When either of the three symmetries, S y , R y , or P is imposed, 
the s.p. Routhian acquires a specific block-diagonal form, and diagonalization of smaller 
matrices results in a faster execution time. In version (v2.08i), this is implemented for 
S y and/or P imposed, but in case R y alone is imposed it is not implemented yet. The 
three T-simplexes and the time reversal cannot be used in a similar way, because they 
are represented by anti-linear operators. However, when the time reversal and simplex 
or signature are imposed, the Kramers degeneracy allows for diagonalization of matrices 
only in one simplex or signature, which reduces the numerical effort by half. The code 
stores the Skyrme local densities on a Cartesian mesh, and takes advantage of imposed 
symmetries, both represented by linear and anti-linear operators, to reduce the number 
of lattice points on which the densities are calculated. 

Arbitrary s.p. configurations can be required in the code HFODD, and the way of 
selecting them is related to the pattern of imposed symmetries. If P and R y are imposed, 
then the s.p. states are divided into four groups characterized by the parity and signature 
eigenvalues. Those groups are referred to as parity- signature blocks. The code does not 
put the particles just on the energetically lowest s.p. states, but the user must specify the 
numbers of particles in each block separately. Then, in each block, the lowest levels are 
occupied, which gives the reference configuration, from which the particle-hole excitation 
are counted afterwards. The same is implemented for imposing either of P, R y , or S y 
alone, only that in these cases there are two blocks instead of four. They are called parity, 
signature, and simplex blocks, respectively. If no symmetries are imposed, obviously 
there is only one block. The three T-simplexes and the time reversal do not play any 
role in this mechanism, because anti-linear operators do not provide quantum numbers. 
In the following, the parity-signature reference configuration is specified as {n ++ , n + _, 

n_ + , n ), where, e.g., n + _ is the number of occupied states of parity +1 and signature 

—i, etc. Similarly, the parity, signature, and simplex reference configurations are given 
as {n + , States in parity-signature blocks are denoted as, e.g., 18-+, which refers 

to the eighteenth state of negative parity and positive signature, counting from one in 
energetical order. Sometimes, a pair of states, like 18-+ and 18 — , is denoted as 18-+. 
The notation for the parity, signature, and simplex blocks is similar, e.g., 28-. Particle- 
hole excitations are defined by selecting from which state in which symmetry block a 
particle is removed and onto which state in which block it is put. This is denoted as, 

e.g., (18- H — ►19-+) for the parity-signature case and, e.g., (28 >29-) for the parity, 

signature, or simplex imposed. The use of symmetry blocks makes it easier to follow a 
given configuration diabatically, because crossings of levels belonging to different blocks 



including the trivial subgroups. 



3.6. THE PROGRAM hfodd (V2.05C) 



41 



Table 3.3: Complete list of imposed (I) or non-imposed (N) symmetries allowed in the 
code hfodd version (v2.08i). The first column enumerates all the possible symmetry 
patterns, while the second column enumerates non-equivalent subgroups (see text). 
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do not cause any confusion as to which of them should be occupied. If two levels from the 
same block are about to cross near the Fermi surface, the user may require that a state 
with a lower/higher value of a certain s.p. observable be occupied. In the present context, 
such diabatic blocking concerns primarily the projection of the s.p. angular momentum 
on the cranking frequency vector. 

The code allows for TAC calculations, and there are no restrictions on the direction 
of the angular frequency vector, except for those following from the imposed symmetries. 
According to the Kerman-Onishi theorem |Ker81aj . in TAC calculations all self-consistent 
solutions have their spin and cranking frequency vectors, J and u>, parallel. In practice, 
J approaches u) very slowly in the course of the iterations, and up to several thousands 
of iterations are necessary to obtain a converged result. To avoid this inconvenience, a 
quadratic constraint on the vector product, u> x J, has been introduced, constraining that 
quantity to zero. However, this does not speed up the convergence significantly, and may 
easily lead to divergencies. In version (v2.08i), a different prescription was implemented, 
which is much more stable and substantially reduces the required number of iterations. 
Namely, in each iteration, the rotational frequency vector is explicitly reset to be parallel 
to the current total angular momentum vector, while its length, u, is kept unchanged. 
Since uj is not constant in this method, it does not correspond to minimization of any a 
priori given Routhian. However, it is the Routhian with the last direction of uj that will 
have been minimized once the converged solution is found. This method was invented 
after the completion of the present calculations, and was not yet implemented in version 
(v2.05c). This is the only difference between the two versions that matters for the present 
study. 

Some quantities, like the multipole moments, only carry a clear information if cal- 
culated with respect to the center of mass. As far as the symmetry patterns used in 
the code (see Table 13 .3|) are concerned, if neither the parity nor T-parity is imposed, 
the center-of-mass of the solution may move away from the origin of the program frame. 
The code then finds its position and recalculates the quantities in question in the center- 
of-mass frame. It is also useful to compute the spin vectors and all types of multipole 
moments in the intrinsic frame of the nucleus, defined by the principal axes of the electric 
quadrupole tensor. Unless any two of the simplex symmetries, S x , S y , S z , S'J, Sy, Sj, 
are imposed, the principal axes may become tilted with respect to the Cartesian axes of 
the program frame; see Section 13.31 The code then finds the directions of the principal 
axes by diagonalization of the quadrupole tensor, and transforms all necessary quantities 
to the intrinsic system of reference. If the parity is not imposed, the solution may acquire 
non-zero linear momentum. The code computes its value and corrects the mean angular 
momenta for that translational motion. 



Chapter 4 

Properties of the h^j2 valence 
nucleons 

Behavior of the shears and chiral bands is, to a big extent, determined by the rotational 
properties of the valence high-j particles and holes. In the cases studied here, we are 
primarily interested in properties of the /in/2 orbital. In this Chapter, some general 
properties of the hu/2 states are analyzed, first within self-consistent one-dimensional 
cranking method, and then analytically. It is demonstrated that those properties are 
affected by the triaxial deformation and that the concerned valence nucleons are expected 
to align their angular momenta rather stiffly along the principal axes of the triaxial shape. 
It is also argued that pairing correlations may soften those alignments. 



4.1 Symmetries of the non-rotating solutions 

In the following two Chapters, Hartree-Fock (HF) shears solutions in 142 Gd and chiral 
solutions in 130 Cs, 132 La, 134 Pr, 136 Pm are presented. Calculations were carried out with 
two Skyrme- force parametrizations, SLy4 [Cha97aJ and SkM* [Bar82aJ, with and without 
the HF time-odd fields included; see Section 13.41 Since there is no physical indication 
for the breaking of parity in either of the shears or chiral bands, this symmetry was im- 
posed in the calculations. For all considered nuclides, symmetry-unrestricted calculations 
(modulo the imposed parity) were performed to check, whether some other symmetries 
are spontaneously broken in the non-rotating states. Whenever the time-odd fields were 
absent (either switched off or self-consistently zero in the ground state of the even-even 
142 Gd), it was found that the mean field possessed the Dj h symmetry. This was inferred 
from comparison of the symmetry-unrestrained results with those obtained with the D\ h 
symmetry imposed 1 , as described in Section EP1 



1 In HFODD, for configurations with unpaired nucleons, the Dj h symmetry can be imposed on the mean 
field by imposing P, S y , Sj, Sj, and switching the time-odd fields off. 
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4.2 Single-particle PAC Routhians for triaxial nuclei 

All the HF solutions considered in the present work are triaxial, and this property deter- 
mines many of their features. It is instructive to examine their rotations about the three 
principal axes first. This amounts to applying Principal- Axis Cranking (PAC), with both 
parity and signature imposed. For non-rotating states with the time-odd terms absent, 
the mean field is anyway symmetric with respect to the whole D\ h group, and imposing 
either of the three signatures does not change it. Although the associated Slater determi- 
nants can be different for different signatures imposed, due to different signature states 
chosen for unpaired particles, this does not affect the mean field; see Section EP1 for de- 
tails. Thus, it is justified to say that three different rotations of the same triaxial object 
are considered. With the time-odd fields present, the symmetry- unrestrained solutions 
violate all the three signatures. Therefore, imposing any of them necessarily leads to a 
substantially different solution, and different for each of the signatures imposed. How- 
ever, the energies, deformations, and s.p. energies are very similar, and the effects from 
the time-odd fields are not so strong. In the code hfodd, it is technically most convenient 
to perform PAC about the y axis 2 . To examine rotation about each given intrinsic axis, 
the non-rotating solutions were first turned in space, by use of appropriate quadrupole 
constraints, so that the required principal axis coincide with the program y axis. 

Sample self-consistent s.p. PAC Routhians for 130 Cs, 132 La, 134 Pr, 136 Pm, and 142 Gd 
are collected in Appendix|Dj The s.p. levels are numbered so that, e.g., 14 — h denotes the 
fourteenth level of negative parity and positive signature; see Section In non-rotating 
solutions, the s.p. levels originating from the spherical hn/2 orbital have numbers 14-± 
to 19-± (each level comprising two states of opposite signatures), for both protons and 
neutrons. This was checked in calculations with constraints on the quadrupole moment 
by the observation that the s.p. energies of the levels 14— ± to 19-± converge to the 
same value with decreasing (3, and no crossings with other negative-parity levels occur 3 
in function of (3. In the spherical case, the /in/2 s.p. levels unambiuously have numbers 
14— ± to 19— ± in all realistic potentials, because hn/2 is an intruder orbital, well separated 
from all other negative-parity states. 

For cranking about the short axis, the lowest sublevels of the hn/2 orbital (14— ±) split 
strongly in function of u, while the highest (19— ±) do not. For the long axis, it is precisely 
the opposite, while for the medium axis only the intermediate levels split. Routhians for 
cranking about the medium axis are given only in Figs. ID. 61 and ID. 71 for the sample case 
of 132 La. In fact, the intermediate hn/2 levels split for cranking about each of the three 
axes, but only weakly. Slopes of the PAC Routhians, e', translate into the s.p. alignments, 
ji, on the i-th axis according to the well-known formula 

de' 

which holds exactly for fixed mean potentials and approximately if the potential changes 
with rotational frequency due to the HF self-consistency or TRS minimization. Therefore, 



2 Although PAC about the x and z axes is also implemented, conservation of the associated simplex 
or signature symmetries does not yield block-diagonal matrices in the basis used Dob04a . 
3 With the exception of the triaxial minima in 134 Pr and 136 Pm; see Section 
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the lowest hn/ 2 sublevels align their angular-momenta vectors on the short axis, and 
the highest on the long one. The corresponding alignments are approximately equal to 
11/2 h. The intermediate states have non-zero alignments on all the three axes. Although 
these results are based on PAC calculations, they are confirmed in the full Tilted-Axis 
Cranking (TAC) results of Chapters El and |B1 For the D^-symmetric mean fields, the 
alignments obtained from one-dimensional cranking acquire a very specific meaning, which 
is discussed in the next Section. 

4.3 D 2~ symmetric mean field 

In all the present HF solutions without the time-odd fields, as well as in virtually all phe- 
nomenological studies of the magnetic and chiral rotation, the mean-field Hamiltonian, h, 
is symmetric under the D\ h group. In such a case, important properties of the s.p. states 
can be inferred from pure symmetry considerations. Eigenstates of a .D^-symmetric po- 
tential exhibit the two-fold Kramers degeneracy, and this Section is focused on properties 
of a single Kramers pair. In fact, to obtain the results of this Section, it is enough to 
assume that h is symmetric with respect to D%, a subgroup of Dj h which consists of the 
time reversal, the three signature, and the three T-signature operations. 

A complete information about the matrix elements of the angular-momentum operator, 
J, between the states of a Kramers pair, (|^), is contained in the real alignment 
vector, «7 M , and the complex decoupling vector, , of the state 

J» = (ji\J\n) , D» = J\fi) . (4.2) 

The alignment and decoupling vectors, J p ', D^, of the state |/2) are related to J M and D^, 
by Eq. (jB.2|) ; see Appendix[Bl For any Hamiltonian, h, invariant under the D\ group, one 
can choose the states of the considered Kramers pair as eigenstates of either of the three 
signatures, Ri, where i = x,y,z. This results in three formally different pairs, 
that correspond to just three different bases in the same two-dimensional eigenspace of 
h. It is shown in Appendix [Bl that there are only three independent real parameters 
that determine all components of the alignment and decoupling vectors associated with 
the three examined pairs. These are the "diagonal" components J£ x , Jy y , Jp , where the 
lower indices refer to the Cartesian components of J. One has 

J"* = (Jt, 0, 0) , = (0, Jf, -Up) , (4.3) 



J»y = (0, J£y, 0) , = (-U£ x ,0, Jp ) , (4.4) 

J"« = (0, 0, J?) , = (jp, -iJf,0) . (4.5) 

Invariance of h under D\ does not restrict the values of Jp, Jy y , Jp. In the often 
considered case of axial symmetry, say with respect to the z axis, 

, (J^,J^,J^ 1 ) for Jj 1 " = 1/2 

{J x % J y , J z ') - <; ( 0j0ji/ /*ii) for = 3/2, 5/2, ... (46) 
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The component J^" is quantized due to the axial symmetry, and the still arbitrary com- 
ponent J^ x is usually called decoupling parameter, see Appendix |B] for details. 

Suppose that the states in the considered Kramers pair have no coupling to other states 
through the angular-momentum operator. Under this assumption, the 3D cranking for this 
pair becomes a two-dimensional diagonalization problem, that can be solved analytically. 
Of course, it concerns the non-selfconsistent cranking, where the Hamiltonian h in the 
Routhian h! is fixed. For a degenerate Kramers pair, h reduces to its eigenvalue, e. Taking, 



e.g., the states \/i z ), 


\fi s ) as 


a basis 


one 


h' = h- 


ujJ = 


" e " 

e 





UJ X J£ X + iUJ y Jy V -U Z J^ 



(4.7) 



For the eigenstates of h', one obtains the following components of the mean angular 



momentum vector, J, 



They differ by sign for the two eigenstates. This dependence has a singularity at u = 0, 
and is not linear in the general case. Consider the following two extreme situations. If 
J£ x = Jy y = J^ z = g, then Jj = guji/ou, and J always orients itself along u>, already for 
infinitesimal u. This can be called soft alignment. On the other hand, when only one of 
the parameters J£ x , Jy y , J^ z is non-zero, say Jj , then Jj = Jj J o~ij, and J is independent 
of u) (unless u = 0). We call this stiff alignment on the j-th axis. 

In axial nuclei, precisely one two-fold degenerate substate of each deformation-split 
j-shell has J^" = 1/2 and J^ x ^ 0, which represents the soft alignment. According 
to Eq. ()4.f)j) . all other necessarily have a vanishing decoupling parameter, and are thus 
rigidly aligned with the symmetry axis. For prolate shapes, the lowest-energy substate has 
J^" = 1/2, and is soft, while for oblate shapes it is the highest substate. In triaxial nuclei, 
values of the parametrs Jf*, where i = s,m,l corresponds to the short, medium, and long 
principal axes, are equal to the alignments obtained from the one-dimensional cranking 
about the three axes. Indeed, for cranking about the axis i, the s.p. states are eigenstates 
of Ri. Taking into account the PAC results of the previous Section, one can see that for 
the lowest hn/ 2 substates of a triaxial nucleus only J^ s is non-zero, while for the highest 
substates only Jf l does not vanish. These alignments are thus stiff. The Routhians of 
these states are not much curved, which confirms that their angular-momentum coupling 
to other states is rather weak. Note that there are no states with stiff alignment on the 
medium axis. The response to rotation of the middle hu/ 2 substates is soft, because 
all the three parameters, J,^" 1 , J^ 1 , are non-zero, and complicated, because these 
states interact with one another, what can be seen from the strong bending of their PAC 
Routhians. 

Although the presented results come from a non-selfconsistent cranking (with fixed h), 
the full HF calculations in Chapters El and El show that at least the deformation does not 
change much with rotational frequency, particularly for the chiral bands. Therefore, the 
non-selfconsistent cranking should be a good approximation here. In both the magnetic 
and chiral bands, the active high-j particles occupy the lowest substates of a given orbital, 
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while the valence holes occupy the highest substates. The principal conclusion of the last 
two Sections is, therefore, that those valence particles and holes align their individual 
angular momenta on the short and long axes of the triaxial shape, respectively, and that 
those alignments are stiff. These results, based mostly on symmetry considerations, agree 
well with the full HF TAC results of Chapters El and According to the presented 
considerations, the triaxial deformation should affect the blades closing in the shears 
bands. Indeed, for prolate or oblate shapes, respectively the lowest particle or the highest 
hole would rather have a soft alignment. 



4.4 Influence of the pairing correlations 

Besides deformation, pairing correlations may affect properties of the shears and chiral 
mean-field solutions. Since pairing is not included in the present HF calculations, this 
Section aims at least at a model study of its possible role. It is well known that pairing 
reduces the collective moments of inertia. Here, the interest is rather in the influence that 
pairing may exert on the alignment properties of the valence particles and holes. 

Non-selfconsistent one-dimensional cranking about the y axis is considered for a system 
of two Kramers pairs, (\/J, y ), \p> y )) and (\v y ), \v y )), whose states have good signature R y . 
Let e M < e u be the energies of those two pairs, taken at uj y = 0, and g y be the matrix 
element of the y-th component, J y , of the angular-momentum operator between \/j, y ) and 
\v y ). It is assumed for simplicity that the considered states have zero alignment on the y 
axis at u y = 0. In the basis of the states \[i y ), \fi y ), \v y ), \v y ), the no-pairing Routhian for 
this model takes the form 



ti 



ei 





-u y g y 








ei 





LO y g y 


u y g y 
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u y g y 





e 2 



(4.9) 



The interest here is in the alignment, j y , of a single particle that at u y = occupies the 
state l/iy). By diagonalization of (j4.9|) and, e.g., by using formula (|4.1|) . one obtains 



Jy 



2u y g 2 y 



+ Auj 2 y gl 



(4.10) 



This dependence is traced in Fig. 14.11 for e M = —1 MeV, e v = +1 MeV and g y — 2 h. 
Inclusion of pairing in the simplest fixed-delta form leads to the quasi-particle Routhian 
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Figure 4.1: Alignment of the canonical particle with v 2 = 1 for sample values of the 
Fermi energy, A, and for the case with no pairing (A = 0), obtained from one-dimensional 
cranking in a simple model (see text). Values of A and A are given in MeV. 

For a qualitative study, a gap parameter of A = 1 MeV is taken. In order to talk about 
an unpaired particle in the Hartree-Fock-Bogolyubov formalism, one has to consider a 
one-quasiparticle excitation and identify the particle in question with the canonical state 
with occupation v 2 = 1. Precisely one such state must appear in the canonical spectrum 
for one-quasiparticle excitation. For u y = 0, that state necessarily coincides with some 
eigenstate of the particle-hole Routhian (j4.9|) . but it is not so for u y ^ 0. In other words, 
in rotating solutions pairing introduces extra admixtures to the blocked particle, that 
can change its alignment properties. To describe a particle initially put on \fi y ), in the 
following such a quasi-particle is blocked that for u y — > the fully occupied canonical 
state continuously pass onto \ii y ). The model was solved numerically. Figure HH] shows the 
alignment, j y , of the considered canonical particle in function of the rotational frequency, 
u y , for a few sample values of the Fermi energy, A. If A <C e M , the result is obviously 
identical to that without pairing. As the Fermi energy increases towards A = 0, the 
process of alignment gets faster. It is illustrated for A = — 1 MeV and A = —0.1 MeV. If 
A = 0, already for infinitesimal rotational frequencies a non-zero value of j y is achieved. 
For this value of A, the Fermi level lies exactly at half way between e M and e u . For A > 
the fully occupied canonical particle initially anti-aligns its angular momentum, and only 
for higher rotational frequencies the positive j y is regained. This is not marked in Fig. 14.11 
because such a situation does not seem to be the case, neither for the magnetic nor for 
the chiral rotation. 

Although no quantitative conclusions can be drawn from this simple model, it becomes 
clear that inclusion of pairing correlations may render the valence particles and holes more 
susceptible to rotational alignment. Together with the reduction of the total moments of 
inertia, it may affect the shears-rotor competition in the magnetic bands and the value of 
the chiral critical frequency, as well as the structure of chiral bands; see Chapter |21 



Chapter 5 

Hartree-Fock shears results in 142 Gd 



In the present work, the Tth\y 2 vh~y 2 single-particle (s.p.) configuration in 142 Gd is taken 
as an example for the study of the shears mechanism. It is particularly simple to handle, 
because it involves only the well-separated hii/% intruder orbital, and can be considered 
a model shears configuration. Additionally, a magnetic band supposed to be built on 
that structure has been observed experimentally in 142 Gd. The presented results show 
the first fully self-consistent Hartree-Fock (HF) Tilted-Axis Cranking (TAC) solutions, 
and confirm the important role of the shears mechanism. However, the agreement with 
experiment is not yet satisfactory. Possible reasons are discussed, but further research is 
needed to clarify this point. A preliminary report on these results was given in [OlbQEl. 



5.1 Previous studies in 142 Gd 

Magnetic bands in 142 Gd were first investigated by Sugawara et al. |Sug97a|, but the 
present analysis is based on the most recent EUROBALL III measurements by Lieder et 
al. |Lie02a| . In that experiment, four magnetic bands, denoted DB1, DB2, DB3, DB4, 
were observed, and configurations of ^h\y 2 vh~^ 2 , irh\y 2 gy 2 ^u/2> 7r ^ii/2 v ^hi/2 anc ^ 
7r ^ii/2fi'f/2 v ^hi/2 were assigned to them, respectively. The present work focuses on the 
DB1 structure. Figure l5~Tl shows a partial level scheme of 142 Gd, including the DB1 and 
the ground-state band. Two Al = 2 bands built on the ^^\\j2 anc ^ v ^hx/2 configurations 
|Sta88aj are also displayed, because if the ith\y 2 and vh~^, 2 excitations are combined 
together, they yield the ith\y 2 vh^, 2 configuration of the band DB1. The ith\y 2 and 
v ^Hxi2 bandheads are 10 + isomeric states, which points to the parallel coupling of the 
angular momenta of the two h\\i 2 proton particles and neutron holes. Such a coupling is 
characteristic for the shears mechanism; see Section 12.11 



It was discussed in Section 12.11 that deformation has a strong influence on the shears 
mechanism. There is no direct experimental data on the deformation of 142 Gd, but several 
theoretical predictions for the ground state are available. The Extended Thomas-Fermi 
method with the Strutinsky Integral |Abo92aj gives (3 2 = —0.21, the Folded Yukawa po- 



tential combined with the Droplet Model [Mol95aJ yields f3 2 = —0.156 and the Relativistic 
Mean Field [Lal96aJ gives (3 2 = —0.158. These theoretical calculations were limited to 
axially-symmetric shapes. In the present work, Total Routhian Surface (TRS) calcu- 
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Figure 5.1: Partial level scheme of 142 Gd. The magnetic dipole band DB1 |Lie02aj is 
shown together with the ground-state sequence and the two AI = 2 bands Sta88aJ built 
on the vr/i 2 1//2 and vh^ 2 configurations. 
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Figure 5.2: Energy surface for the non-rotating ground state in 142 Gd, obtained from the 
TRS calculations taking into account the pairing correlations. 
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lations allowing for triaxial deformations were performed with the pairing correlations 
included; see Appendix O for details. The obtained energy surface is shown in Fig. 15.21 
The minimum corresponds to a triaxial shape with 7 « 30°. At 7 = 0° and 7 = 60° 
there are saddle-points, and the one for 7 = 60° has lower energy, in agreement with 
the previous predictions of oblate shapes. The triaxial minimum is rather shallow in the 
7 direction, so one can talk about 7-softness. Both the minima and the saddle-points 
have the absolute deformation of (3 ~ 0.181, placing the TRS prediction between those 
of |Mol95a[ ILal96aj and that of [Abo92aJ. For the magnetic bands in 142 Gd, phenomeno- 
logical TAC calculations allowing for triaxial shapes were performed in the experimental 
paper |Lie Q2aj . by using the code of Frauendorf; see Section 13.51 Energy minimization 
gave (3 ~ 0.15 and 7 ~ 60° for all the four bands, but (3 ~ 0.10 and 7 = 60° was taken 
for the cranking description of DB1 and DB3. A satisfactory agreement with experiment 
has been obtained for both the energies and the B(M1)/B(E2) ratios. 

5.2 HF solutions in 142 Gd 

This Section describes the HF ground-state and shears solutions in 142 Gd. Since this is 
the first application of the Skyrme-HF TAC method developed in this work, below it 
is explained step-by-step how those solutions were found. Refer to Section 13.61 for the 
treatment of symmetries and s.p. configurations in the code hfodd. 

The HF calculations in 142 Gd were performed with two Skyrme parametrizations, 
SLy4 [Cha9Zij and SkM* iBar82a ; see Appendix O for all details. To investigate the 
role of the HF time-odd densities and fields, three different sets of time-odd fields, N 
(no time-odd fields), G (time-odd fields responsible for a gauge invariance of the force), 
and T (all Skyrme time-odd fields), were included; see Section "3~4l The study of the 
shears mechanism was preceded by a search for the ground-state minimum. Only the 
parity was imposed as an unbroken symmetry. Several parity reference configurations (see 
Section l3~T)|) were considered that seemed reasonable from the standard Nilsson diagrams, 
with no particle-hole excitations required. Cranking constraint was applied in the first 
iteration to make the calculations symmetry-unrestricted, and then released, as described 
in Section 13*31 The lowest minimum was found for the parity configuration listed in the 
first row, right column of Table 15.11 With the SLy4 force, that solution corresponds to a 
HF vacuum, and for SkM* there is one empty positive-parity state under the last occupied 
negative-parity state. It is known that the order of s.p. levels may vary from force to force. 
Deformation of (3 = 0.19, 7 = 40° for SLy4 and (3 = 0.18, 7 = 39° for SkM* was obtained. 
These values of (3 are very close to the TRS result for the ground state in 142 Gd. 

It should be expected that the ground state of the even-even 142 Gd with quadrupole 
deformation is symmetric with respect to the whole D\ h group, not only the parity. 
If that was the case, each Kramers pair of a given parity would comprise two states 
of opposite signatures, and the considered parity configuration would translate into a 
parity-signature configuration given in the first row, left column of Table 15.11 Indeed, 
-D^-imposing calculations with such a configuration were performed, and an identical 
result was obtained; see Section 13.31 

The Z^-imposing solution was then used to examine one- dimensional rotations about 
the three principal axes; see Section 14.21 Figure ID. II shows the s.p. Routhians from the 
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Table 5.1: Parity-signature and parity single-particle configurations (see Section RTTTTf) for 
protons (it) and neutrons (z/) used to obtain the solutions in 142 Gd at zero cranking 
frequency 









Parity- S ignat ure 




Parity 








Ref. Exc. 


Ref. 


Exc. 


142 Gd 


SLy4 


TT 


(18,18,14,14) 


(36,28) 




Ground 


SkM* 


V 


(21,21,18,18) 


(42,36) 




142 Gd 


SLy4 


TT 


(18,18,14,14) (14 >15— ) 


(36,28) 


(28 >29-) 


Shears 


SkM* 


V 


(21,21,18,18) (18- + ->19-+) 


(42,36) 


(36 >37-) 



Principal-Axis Cranking (PAC) calculations with the SLy4 force and no time-odd fields. 
As it was discussed in Section 14.21 the s.p. levels 14-± to 19-± (see Section 13 .6J) can 
be recognized as belonging to the hn/2 orbital for both kinds of nucleons. One can see, 
therefore, that in the considered HF solution there are two hu/ 2 proton particles occupying 
the lowest two-fold degenerate substate of the orbital and two /tn/2 neutron holes on its 
highest substate. 

The considered shears band was calculated with only the parity imposed. In the 
language of the parity configurations, the ^11/2 orbital comprises states from the one 
described by 27- to that described by 38-, both for neutrons and protons. With no time- 
odd fields, the Trh^^ vh~^, 2 shears configuration was obtained from the ground state 
by exciting one proton from the 28- onto the 29- state, and simultaneously promoting 
one neutron hole 1 from 37- onto 36-, see the second row of Table l5~Tl This led to a 
converged solution for zero cranking frequency. Still with the time-odd fields switched 
off, a similar test as for the ground state was performed, showing that the mean field in 
the TTh^^ vh~y 2 solution breaks none of the D% h symmetries with the parity-signature 
configuration as specified in the second row, left column of Table 15.11 

To obtain the two-dimensional TAC solutions, a cranking frequency of 0.05MeV/fi 
with equal components on the short and long axes of the nucleus was applied, starting 
the iterations from the previously converged non-rotating state. For the sake of sim- 
plicity, the mean field obtained with the Dj h symmetry imposed was used as a starting 
point, which guaranteed that the principal axes initially coincided with the program axes; 
see Section 13.31 As expected, a perpendicular coupling of the valence-nucleon angular 
momenta was obtained. The band was calculated with the cranking-frequency step of 
0.05 MeV/h, starting the iterations for each subsequent value from the converged solution 
for the previous one. 

The two occupied proton hn/2 states are marked with filled circles in Fig. ID.2l showing 
the s.p. Routhians for the shears band with N time-odd fields. These are just the lowest 
two /in/2 states with positive alignments on the cranking vector, as reflected by their 
negative slopes. At low frequencies, the two states have numbers 27- and 29-, as discussed 
previously. At uj « 0.15 MeV/h, there is a crossing of the down-sloping 29- Routhian 



1 Since the code hfodd works in the language of particles rather than holes, this was actually done by 
removing one neutron particle from 36 -, and putting it onto 37 -. 
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with the up-sloping 28- Routhian. For the configuration to be followed diabatically, the 
proton particle has to be moved from 29- onto 28-, and one deals from now on with a 
vacuum state in the negative-parity block. The analogous holds for the two neutron holes. 

Obtaining the HF TAC solutions in the presence of the time-odd terms of the Skyrme 
mean field is more cumbersome than without those terms. In general, the inclusion of the 
time-odd fields splits the Kramers pairs even in non-rotating states, which may change 
the order of the s.p. levels. In particular, some empty states from above the Fermi surface 
may come below the highest occupied states, which formally corresponds to particle-hole 
excitations. Since one cannot foresee which such excitations are a priori needed, it is 
difficult to find the s.p. configuration corresponding to the required physical state at zero 
cranking frequency. It is easier to proceed in case the signature or simplex symmetry is 
imposed, because such symmetry provides an additional quantum number, which allows 
to distinguish between s.p. states. Moreover, practice shows that calculations with the 
time-odd fields are prone to divergencies at low cranking frequencies. 

In order to circumvent the above problems, calculations with the time-odd fields were 
not started from the zero cranking frequency. Instead, iterations for the G and T fields 
were restarted at each cranking frequency from the corresponding converged solution with 
the N fields, with the s.p. configuration unchanged. Because of shifts in the s.p. energies 
introduced by the time-odd fields, this did not work in cases when levels cross either 
due to those fields, or in function of the rotational frequency. Yet, converged results 
were obtained for higher frequencies, where the proton and neutron configurations are 
already the s.p. vacua in negative parity. These fragments of bands were then extended 
towards lower frequencies step by step, by changing the configuration by hand whenever 
a crossing was about to occur, or by using the diabatic blocking on the total alignment on 
the cranking vector; see Section |3~H1 In spite of that, convergence could not be achieved 
in some cases. See Fig. 15.41 for the frequency ranges, in which the G and T solutions did 
converge. 

5.3 The shears mechanism in 142 Gd 

The key point of the present study is to examine the shears mechanism in 142 Gd. It is 
particularly interesting to see how the aligning angular momenta of the valence nucleons 
and the collective rotation compete in generating the total spin. These issues are addressed 
in this Section. 

Below, spins of the two valence neutron hu/2 holes are defined as the spins of the two 
empty hu/2 states, taken with the minus sign. To ensure that the total neutron spin is 
equal to the sum of the two holes' and the core's contributions, the core is defined as all 
the occupied neutron states plus the two empty h u /2 states. The division of the proton 
s.p. states into the two valence h u /2 particles and the remaining core is obvious. 

The angular momentum vectors of the valence nucleons and of the proton and neutron 
cores are shown in Fig. 15.31 for several values of the rotational frequency. The angular 
momenta of the two hu/2 proton particles are nearly parallel (stretched coupling), and 
their vector sum (particle blade) has length of approximately 10h. The same holds for 
the neutron holes. The stretched coupling and the length of the two blades is retained up 
to high cranking frequencies. In the non- rotating state (bandhead), the particle and the 
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Figure 5.3: Angular momentum vectors of the valence nucleons and of the proton and 
neutron cores along the shears Tr/i 2 ^ , 2 vh^ 2 band in 142 Gd. HF results with the SLy4 and 
SkM* forces are shown for the N, G, T time-odd fields included. For the proton particles, 
the arrow anchored at the origin stands for the energetically lower of the two concerned 
s.p. states, while the inverse holds for the neutron holes. Arc of radius 10h is drawn in 
each panel, and the angle between the proton and neutron blade is given. Rotational 
frequencies are given in MeV/h. 
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Figure 5.4: The (3 and 7 deformations along the shears vr/i 2 -^ vh xi i 2 band in 142 Gd. HF 
results with the SLy4 and SkM* forces are shown for the N, G, T time-odd fields included. 



hole blade aligns with the short and the long intrinsic axis, respectively. In the absence 
of collective rotation, that perpendicular orientation yields the total bandhead spin of 
about 14h. With increasing rotational frequency, the two blades align towards each other 
(shears closing). These results are perfectly consistent with the shears scenario proposed 
in the literature and discussed in Section ETT1 At uj = 0.75 MeV/h, the angle between the 
blades is equal to 53° for SLy4 and 50° for SkM*. At that point, the angular momentum 
gain due to the shears closing (with respect to the bandhead spin) reaches about 5h. Note 
that the maximum gain, corresponding to the complete closure, would be not more than 
6h for two blades of 10^, each, initially forming a right angle. 

In competition with the shears mechanism, also the core generates its angular mo- 
mentum. Both the proton and neutron cores tend to align their spins with the long axis, 
because the moment of inertia with respect to that axis is larger than with respect to 
the short axis. This is expected from the irrotational-flow formula (|2.9|) for 7 > 30°, and 
can be also checked by considering one-dimensional cranking about the two axes, like in 
Section 16.31 Interestingly, the contribution from the neutron core is a few times smaller 
than that from the proton core. At uj = 0.75 MeV/h, the spin generated by the two cores 
equals ~ lOh, two times more than the contribution from the shears mechanism. This 
is in accord with the estimate of Macchiavelli and Clark |Mac99aj (see Section 12. 1|) that 
deformations smaller than e w 0.12 are needed for the shears mechanism to dominate. 
It can be seen from Fig. 15.31 that the behavior of the shears blades and of the two cores 
is very similar for the SLy4 and SkM* Skyrme parametrizations, and independent of the 
time-odd fields included. 
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At io = 0, the shears solution has the deformation of (5 ~ 0.17, 7 ~ 45°, with negligible 
differences between the forces used. Thus, the nh^^ vh xl , 2 excitation slightly reduces 
the value of (3 as compared to the ground state. The deformation further diminishes 
with increasing rotational frequency, which is illustrated in Fig. 15. 4| and falls below 0.15 
with uj approaching 0.75 MeV//L Such a decrease of deformation is characteristic for the 
shears bands; see Section ITT1 In the present results, however, the solution simultaneously 
becomes more and more triaxial, as it can be seen from the evolution of 7 shown in the 
same Figure. The closing blades composed of the two pairs of high-j nucleons have quite 
a strong polarization effect on the core. 



5.4 Comparison with experiment 
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Figure 5.5: Calculated energies for the Tth\y 2 vh^ X j 2 shears band in 142 Gd (lines), com- 
pared with the experimental data for the DB1 band (points). Hartree-Fock results with 
the SLy4 and SkM* forces are shown for the N, G, T time-odd fields included. The energy 
scale refers to the experimental energy measured with respect to the ground state. For 
the TV and G fields, the calculated curves are offset in energy so that they coincide with 
the experimental bandhead at spin 16^. For the T fields, the offset is arbitrary. 



In Fig. 15.51 the energies calculated for the Tth xl ^ 2 vh~y 2 band are compared with the 
experimental data for the band DB1. Clearly, in the calculations the spin is generated at 
a too low energy expense in the vicinity of the experimental bandhead. Since the shears 
closing is rather slow in the HF solutions, the discrepancy should rather be attributed to 
an overestimated collective inertia. The effect is most pronounced for the G set of the 
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time-odd fields, which systematically yields the largest values of the moments of inertia; 
also for the N = 75 isotones studied in Section 16.51 Apart from the time-odd fields, 
the moment of inertia depends strongly on the deformation and on the intensity of the 
pairing correlations; see Section 14.41 Since the nh^^ vh^y 2 shears configuration is a 
two-quasiparticle state in each kind of nucleons, the pairing effects should be significantly 
weakened, and one would expect the pure HF approximation to provide a correct de- 
scription. From the theoretical point of view, there is no apparent reason to surmise that 
the calculated deformation is wrong, all the more since the predictions for the ground 
state are consistent with the TRS results. It must be conceded, however, that the quoted 
estimates of Mol95a, Lal96aJ give a considerably lower ft for the ground state. 

Whether the combined effects of the lack of pairing in the present calculations and 
of a possibly overestimated deformation may account for the observed deviation from 
the experimental data should be a subject of further research. Certainly, there are still 
several unclear points about the theoretical description of the magnetic bands in 142 Gd, 
particularly concerning the deformation. In |Lie02aj . a considerably smaller e had to be 
used to reproduce the DB1 and DB3 bands than it came from the energy minimization. 
From the level scheme in Fig. 15. H it is clear that the bands built on the 7r/i 2 1//2 and vh^ 2 
configurations are rather well-deformed structures, similar to the ground-state band. Why 
would the combination of these two excitations diminish the deformation so much? It 
should be noted here that recently Pasternak et al. |Pas04aj obtained a surprisingly 



good description of the magnetic bands in 142 Gd within a very simple classical model. It 
would be instructive to see whether such solutions have their self-consistent mean-field 
counterparts. 
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Chapter 6 

Hartree-Fock chiral results 



Investigation of the chiral rotation in this thesis was directly stimulated by the first 
experimental observations of candidate chiral bands, in 130 Cs, 132 La, 134 Pr, and 136 Pm 
[StaO^j. These nuclei were naturally chosen for the analysis in this work. An exploratory 
study performed within the Principal- Axis Cranking (PAC) revealed that a very appealing 
and accurate model for the chiral rotation can be given in the classical framework of two 
gyroscopes coupled to a triaxial rotor. That model is presented in Section lowl and leads 
to one of the principal conclusions of this dissertation, that chiral rotation can only exist 
above some critical value of the rotational frequency or spin. This result is confirmed 
by the full Skyrme-Hartree-Fock (Skyrme-HF) solutions that are presented later in this 
Chapter, and constitute the first fully self-consistent proof of the chiral symmetry breaking 
in rotating nuclei. However, the chiral bands could be obtained only in 132 La from among 
the studied isotones. A brief report on the results obtained in 132 La is given in jQTb04a , 
which is a revised version of the e-print published earlier in Olb02bJ. 



6.1 Previous studies in 130 Cs, 132 La, 134 Pr, and 136 Pm 

As it was discussed in Section \2.2\ the Tch^^ vh^y 2 configuration is in the literature 
commonly adopted for the chiral doublet bands in the N = 75 isotones. Level schemes 
of those bands are displayed in Figs. 16.11 and 16.21 Experimental data were taken from 
the most recent measurements, namely from |KoiQ3al for 130 Cs, from |Sta02at ITim03a( 



ICro04a| for 132 La, from jStaOlbl |Rob03a for 134 Pr, and from |Har01a| in the case of 
i36p m rj^g names G f bands in 132 La, Bl and B3, are reproduced from [Gro04aj. Spin 
assignments are still tentative, apart from the two lowest bands in 132 La. Note that the 
bands referred to as yrast in the present work are only yrast in the positive parity. The low- 
spin regions of those bands are difficult to explore using 7-spectroscopy, because transition 
energies become very low. However, recent highbrow experiments in 132 La |Tim03a| and 



134 Pr |Rob03a| suggested that the bandhead spins are 7 + , which is consistent with the 
orthogonal angular momentum coupling of the hu/2 proton particle and neutron hole. 
See Tables 12.11 and 12.21 for a review of experimental and theoretical investigations of the 
chiral rotation in the N = 75 isotones. 

For all of nuclei 130 Cs, 132 La, 134 Pr, and 136 Pm phenomenological Tilted- Axis Cranking 
(TAC) calculations were performed in [StaOlbj . and in [DimOOa also for 134 Pr and in 
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[HecOlaj for 136 Pm. All these calculations were performed by using the TAC code of 
Frauendorf. The e deformations of 0.16, 0.175, 0.175, 0.195 and triaxialities 7 of 39°, 
32°, 27°, « 26° were found, respectively for 130 Cs, 132 La, 134 Pr, 136 Pm. Chiral solutions 
were obtained in a limited range of angular frequency. The lower limits correspond to 
the critical frequency, u crit , discussed in this Chapter. Only IDimO Oa^ and [HecOlaj quote 
their values, which are uj cr u = 0.3 MeV/ft for 134 Pr and uj crit < 0.2MeV/h for 136 Pm. 
Comparison with experimental energies is given only for 134 Pr [Sta Olbt IDimOOaj , where 
the mean trend is reproduced. 

For all the considered N = 75 isotones, calculations within the Particle-Rotor Model 
(PRM) were performed as well |Pen03aj . A different version of the model was used 
in [HarOlaj for 136 Pm. The same deformations were used as in the TAC calculations 
described in the previous Paragraph. A satisfactory agreement with experimental data 
was obtained. Probably the most sophisticated implementation of the PRM was employed 
in 132 La jStaOlal ISta02aj and 134 Pr |Sta02aj . Deformations of (3 = 0.23, 7 = 21° and 
(3 = 0.25, 7 = 35° were adopted, respectively, and good agreement was obtained both for 
energies and B(M1)/B(E2). Note that the values of (3 used in the latter calculations are 
significantly larger than in the former ones. 

Recently, a lifetime measurement in 132 La [Gro04aJ provided the first data on the 
absolute values of B (Ml) and B(E2) in the proposed chiral bands. That experiment 
revealed that the band Bl, so far taken for the chiral partner, has the intraband B(E2) 
values an order of magnitude lower than the yrast band and than those predicted by the 
PRM. Moreover, a new band, called B3, was discovered, that cannot be ruled out as the 
chiral partner, either; see Fig. 16.151 Its electromagnetic properties are more similar to 
those of the yrast band. Energetically, it is more distant form the yrast band than Bl, 
but the spin assignments are tentative, and there are some indications |Gro04bj that they 
should be shifted up by one unit. In such a case, the band B3 would become almost 
exactly degenerate with the yrast sequence at spin 18 h. Certainly, further experimental 
research is needed to clarify this point. 



6.2 HF minima in the N = 75 isotones 

The present HF study of 130 Cs, 132 La, 134 Pr, and 136 Pm is limited to the chiral ivh\y 2 vh~i/ 2 
configuration. This Section describes the structure of the minima that were found for that 
configuration before application of cranking. 

Calculations were performed with two Skyrme parametrizations, SLy4 [Cha97aj and 
SkM* [Bar82jiJ; see Appendix O for all details. To investigate the role of the HF time-odd 
densities and fields, three different sets of time-odd fields, N (no time-odd fields), G (time- 
odd fields responsible for a gauge invariance of the force), and T (all Skyrme time-odd 
fields), were included; see Section I3~4l The nh^^ vK^^ configuration is energetically 
most favored in the mean field for the single hx\ii proton particle occupying the lowest 
substate of the orbital, and the neutron hole occupying the highest substate. Only such 
a case was considered, and the HF solutions were sought much like those in 142 Gd, see 
Section l5~2| with the parity imposed. It was inferred from comparison of the symmetry- 
unrestricted and .D^-imposing calculations that the D\ h symmetry was not spontaneously 
broken in all solutions where no time-odd fields were included. Table l6~T1 gives the parity 
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Figure 6.1: Partial level schemes of 130 Cs |Koi03a| . 134 Pr jStaOlbi IRob03a| and 136 Pm 
|Har01aj . showing the yrast positive-parity bands (to the right) and supposed chiral part- 
ners (to the left). 
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Figure 6.2: Partial level scheme of 132 La, showing the yrast band (middle) and its two 
candidate chiral partners. The band Bl is the previously known partner |Sta02aj . while 
B3 is the recently discovered third band |Gro04aj . 

and parity-signature configurations (see Section l3~K|) that were found for the studied states. 

In 130 Cs and 132 La, deformations of (3 = 0.24 — 0.26 and 7 45° — 49° were obtained, 
depending on the force and time-odd fields. In the N = 75 isotones, 25 neutrons are left 
above the N = 50 spherical shell. Within the considered configuration, 11 of them are put 
on the /in/2 orbital, and the remaining 14 occupy the positive-parity single-particle (s.p.) 
states above iV = 50. Presumably, they entirely fill the g 7 / 2 and d 5 / 2 orbitals, and thus 
exert no shape-driving force, although the positive-parity orbitals are obviously mixed by 
deformation. 

In the Z = 55 nucleus 130 Cs, one of the 5 protons outside the Z = 50 shell occupies 
the /in/2 orbital, and the other 4 the #7/2 orbital. Addition of two protons ( 132 La) also 
results in their location on that orbital. Although there obviously exist some mixing 
between the #7/2 and d 5 / 2 proton orbitals, they seem to be quite well separated and can 
be easily recognized among the s.p. Routhians, see Figs. ID. 51 and ID.6[ ID. 71 thanks to 
their alignment properties on the short and long axes. The $7/2 orbital is composed of the 
levels 13+±, 14+±, 15+±, and 16+± (see Section l3~Bj) . and d 5 / 2 comprises levels 17+±, 
18+±, and 19+±. 

In both 134 Pr and 136 Pm, two minima with the same fth\y 2 vh~^ 2 configuration were 
found, which differ by the occupation of positive-parity states. The energetically lower 
minimum has similar positive-parity s.p. structure as in 130 Cs and 132 La. Namely, the two 
protons added to form 134 Pr take the last free substates of #7/2, while the two protons 
added to form 136 Pm occupy d 5 / 2 . This can be traced in Figs. ID.8I and ID. 101 However, 
such solutions correspond to almost oblate shapes of 7 = 53° — 58°. The other minima 
have 7 = 19° — 22°, which is triaxial, but closer to the prolate shapes 1 . It seems that the 
occupied positive-parity neutron g 7 / 2 and d^/ 2 orbitals are now strongly mixed with the 
S1/2 and g?3/2 states from the same magic shell, which may also be the case for protons. For 

1 In the following, those two minima in 134 Pr and 136 Pm are referred to as oblate and triaxial. 
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such deformations, some substates of the negative-parity fi/ 2 or hg/ 2 orbitals from above 
N = 82 come so low that they slip under the highest h\\/ 2 substate, where the valence 
neutron hole is put 2 . The highest hn/ 2 state can be distinguished thanks to its alignments 
on the short and long axes, which are exactly opposite to those of the down-coming lowest 
substates of f 7 / 2 or h 9 / 2 . In Figs. ID.91 and HXTTj the highest hiy 2 substate has numbers 
20-±. Exact deformations in the obtained minima are displayed in Table 16.21 

Table 6.1: Parity-signature and parity single-particle configurations (see Section l3~^j) for 
protons (it) and neutrons (z/) used to obtain the solutions in the iV = 75 isotones at zero 
cranking frequency. 
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6.3 PAC calculations 

Before going to the full HF TAC calculations in the N = 75 isotones, it is instructive 
to have a closer look at the PAC results, already discussed in Section l4~2l They contain 
an important information on both the valence nucleons and on the collective core. The 
results in this Section lead to the formulation of the classical model, that will serve as a 
guideline for a consecutive HF study of the chiral rotation. 

Below, the study of the valence nucleons is carried out for the PAC calculations in 
132 La, performed with the SLy4 force and the G time-odd fields. The conclusions derived 
in this sample case are supported for other nuclei and forces by the full TAC results; 
see Sections 16.61 and 16.81 Figure 16.31 shows the s.p. alignments of the valence particle 
and hole, obtained from one-dimensional cranking about the three principal axes. The 
plot confirms that, at uj = 0, the particle and hole do indeed orient their spins, j p and 

2 There is only one such intruder with the SkM* force, but two with SLy4. This makes the diabatic 
cranking calculations with the SLy4 force difficult due to multiple level crossings. In the following, only 
the SkM* results are shown for the triaxial minima, and only the SLy4 results for the oblate ones. 
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Figure 6.3: Single-particle alignments of the valence hn/ 2 proton particle and neutron 
hole in 132 La, obtained from one-dimensional cranking about the short, medium and long 
axis. The Figure shows the Hartree-Fock results for the SLy4 force and the G time-odd 
fields included. 

j , on the short and long axis, respectively. The response of j p and j to rotation is 
gradual and rather weak, meaning that the s.p. wave functions are strongly constrained 
by deformation (deformation-alignment). These results can be summarized as 

f ~ s s i s + 5J p uj , j h ~ s,i, + 8J h u , (6.1) 

where i s and i\ denote the unit vectors along the short and long axis, respectively. There- 
fore, to a reasonable approximation, the odd particle and hole can be treated like gy- 
roscopes of spins s s and Sj, rigidly fixed along the short and long axes, while the small 
coefficients bj v and 5J h can be incorporated into the total collective inertia tensor, J . 

The same PAC calculations provide the total alignments, J s , J m , Ji, that are plotted 
in Fig. 16.41 for all the considered nuclei and forces. The upbend in the J s curve for 130 Cs 
is caused by the change of slope of the last occupied positive-parity proton s.p. state due 
to its interaction with the first empty one. These are the levels 14++ and 15++ in the 
top-left panel of Fig. ID.5| respectively In 132 La, the 14++ and 15++ states are both 
occupied, and no upbend occurs. Similarly, the upbend in J/ for the triaxial minimum 
in 136 Pm is because of the interaction between the neutron hole state 18++ and the 
particle state 20++, see the bottom-right panel in Fig. ID. Ill The upbend in J s for the 
nearly oblate solution in 136 Pm originates from the smooth crossing of the 17++ and 
18++ proton states, as shown in the top-left part of Fig. ID. 101 This crossing has no 
effect in 134 Pr, because there are two protons less there, and the concerned states are 
both empty. In the nearly oblate minima in 134 Pr and 136 Pm, the medium axis is not well 
distinguished from the long one, and an attempt to crank around it immediately leads to 
solutions associated with the long axis. 

One can see that all the observed bends in the total alignment plots are due to smooth 
level crossings near the Fermi surface. Otherwise, the Ji{oJi) dependence is, to a good 
approximation, linear, like for the rigid rotation. The corresponding slopes give the 
collective total moments of inertia, J S) J m , Ji, with respect to the short, medium, and 
long axis, respectively. They already contain the small contributions, bj v and SJ~ h , from 
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Figure 6.4: Total alignments obtained from one-dimensional cranking about the short, 
medium and long axis in 130 Cs, 132 La, 134 Pr, and 136 Pm. Hartree-Fock results with the 
SLy4 and SkM* forces are shown for the N, G, T time-odd fields included. 
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Figure 6.5: The classical model of chiral rotation. Two gyroscopes of spins s s and si 
rigidly fixed along the short and long axes of a triaxial rigid body characterized by three 
moments of inertia, J s , J mi Ji, associated with the short (s), medium (m) and long (7) 
axis. 

the valence particle and hole 3 . At zero frequency, cranking around the medium axis gives 
a vanishing angular momentum, while those around the other two axes give non-zero 
values, equal to the initial alignments, s s and si, of the odd particle and hole. 

The microscopic PAC results presented so far suggest that the considered system can 
be modeled by two gyroscopes of spins s s and si rigidly fixed along the short and long axes 
of a triaxial rigid rotor characterized by the inertia tensor, J ', whose diagonal components 
are equal to J s , J m) J\. It is instructive to solve the associated problem of motion in the 
classical framework, which is done in the next Section. It was mentioned in Section E~T1 
that for no time-odd fields included in the present HF calculations, the mean potentials for 
the non-rotating states have the D\ h symmetry, which implies that the inertia tensor must 
be diagonal in the principal-axis frame of the mass distribution. With the time-odd fields 
included, perturbative TAC tests were performed, in which to the self-consistent non- 
rotating solution cranking frequency was applied in different directions in order to check 
the response of the mean angular momentum vector. It was found that the off-diagonal 
components of J are negligible. Therefore, the inertia tensor is taken diagonal. 



6.4 Classical model and the critical frequency 

It was argued in the previous Section that a simple model of two gyroscopes coupled to a 
triaxial rigid body is appropriate for the description of rotation in the considered nuclei. 
Below, that model is solved in the classical framework. It possesses several analogies to 
the HF TAC method in its foundations. Yet first of all, it leads to a crucial conclusion 
that chiral rotation can only exist above some critical value of the rotational frequency 
or spin. A simple analytical formula allows for estimating the critical frequency on the 
basis of the standard PAC results. 

3 Since the coefficients Ji, i — s,m,l, are slopes of the Ji(wi) dependence, they should be thought 
of more as the second, or dynamical moments of inertia. In classical mechanics, they would have the 
meaning of components of the inertia tensor. 
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The model system is sketched in Fig. 16.51 A triaxial rigid body models the collective 
core. Its inertia tensor, J , is diagonal with respect to the short, medium and long axes, 
and the corresponding components, J sy J my J7i, also include the small contributions from 
the valence particle and hole; see Eq. (j6.1jl . The particle and hole are represented by two 
ideal gyroscopes of spins s s and si, rigidly fixed along the short and long axis of the core, 
respectively. 

The notion of the ideal gyroscope requires explanation. Consider an axially symmetric 
rigid body with moments of inertia, J\\ and J7j_, with respect to the symmetry axis and 
an axis perpendicular to it, respectively. If such a body spins at angular frequency Q 
around the symmetry axis, its kinetic energy is T = J'\\Q 2 /2. Simultaneously, if the 
body rotates with arbitrarily oriented frequency, uj, the kinetic energy becomes T w = 
( J ± oo 2 sin 2 -Q + J\\{ Locosd + fl) 2 )/2, where $ is the angle between u> and the symmetry 
axis. The ideal gyroscope is the limiting case when J\\ and J7j_ go to zero and Q increases 
so that QJ'w remain constant and equal to a given value of s. Obviously the total spin is 
then equal to s and points along the symmetry axis, while T u becomes infinite. However, 
the difference T u — T = (J7j_ sin 2 $ + J\\ cos 2 ?9)f2 2 /2 + J\\VLuj cost? has a finite limit, 

T = cjscos?? = ws . (6.2) 

This difference can be considered as the kinetic energy of the ideal gyroscope measured 
with respect to the case with no rotation, that is with u> = 0. 

Coming back to the discussed classical model, the total angular momentum, J, of the 
system reads 

J = Juj + s , (6.3) 

where s is the vector sum of spins of the the two gyroscopes. In the absence of potential 
interactions, the Lagrangian of the system is equal to the total kinetic energy, and is given 
by the formula 

L = E kin = -ujJu + ujs . (6.4) 

Taking the laboratory components of u> as generalized velocities, it is easy to check that 
the generalized momenta are equal to the laboratory components of J. This fact allows 
us to write the Legendre transformation [Gol53a j and to obtain the Hamiltonian of the 

system, 

1 R 2 R 2 R 2 

H = uJ-L = -uju> = + ^f- + ^= , (6.5) 

2 2j s 2 l y m 2^7; 

where R = J — s is the angular momentum of the rigid body alone. Note that H is 
analogous to the rotational part (|2.8J) of the PRM Hamiltonian. Since the Lagrangian 
(16.4)1 does not depend explicite on time, Hamiltonian ()6.5)) is a constant of motion, and 
is identified with the total energy, E, of the system. Consider now a particular type of 
the Routh function |Gol53aj (Routhian), H', namely such that no variables undergo the 
Legendre transformation. In such a case, H' = —L, and by rewriting the Routhian in 
terms of the Hamiltonian one obtains 

H' = H-ojJ. (6.6) 



This quantity is similar to Routhian (j3.1())) appearing in the quantum cranking model. 
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Equations of motion for the model can be derived in the following way. The time 
derivatives, d t and &f, of any vector J, taken in the laboratory frame and in a frame 
rotating at angular frequency u>, respectively, are related by the formula d t J = d^J+ujxJ 
[Gol53aJ. But dt J = 0, if J represents the angular momentum, which is conserved in the 
laboratory frame. Thus, one obtains the Euler equations |Gol53aj for the time evolution 
of the angular-momentum vector in the body-fixed frame, 

&tJ = -u> x J . (6.7) 



Thouless and Valatin |Tho62aj showed that analogous equations govern the evolution of 



the mean angular momentum in the Time-Dependent HF theory. 

Apart from the equations of motion, a convenient way to analyze rotation of the 
considered system is through phase portraits, which are intersections of the angular- 
momentum sphere with the energy ellipsoid. Dimitrov et al. carried out such an analysis 
|Dim02aj for a slightly simplified version of the model. Separatrices defining the planar 
and chiral rotational regimes as well as the chiral vibrational regime (see Section 12.21) 
are clearly visible in this approach. However, the mean-field cranking approximation can 
only account for the so-called uniform rotations, in which the mean angular-momentum 
vector is constant in the intrinsic frame of the nucleus, dfJ = 0. Because of that, also the 
present study within the classical model is limited to such uniform motions. The Euler 
equations ()6.7j) now take the form w x J = 0, and require that & and J be parallel. The 
same condition holds for the HF solutions and is known as the Kerman-Onishi theorem; 
see Section 13.31 and |Ker81aj . These equations can be easily solved for the considered 



model. However, to show further analogies with the HF method, here the variational 
principle is employed to find uniform solutions. 

According to the Hamilton's principle, motion of a mechanical system can be found 
by making the action integral, J L dt, stationary. Here, real uniform rotations are sought, 
and all of them obviously belong to a wider class of trial motions with u) constant in the 
intrinsic frame. Within this class, Lagrangian ()6.4)1 does not depend on time. Therefore, 
extremizing the action for the given u reduces to finding extrema of the Lagrangian as a 
function of the intrinsic-frame components of u>. Since H' = —L, the Routhian ()6.6|) can 
be equally used for this purpose. This serves us as a bridge between the classical model 
and quantum cranking theory, where an analogous Routhian is minimized in the space of 
the trial wavef unctions. 

Extrema of H' with respect to the intrinsic-frame components of a; at a given length 
of uj can be found by using a Lagrange multiplier, //, for uo 2 . Setting to zero the derivatives 
of the quantity 

H' + ^fiu 2 = - [(/i - J s )u 2 + J m )uj 2 m + Jfiuf] - {u s s s + msi) (6.8) 

with respect to u s , co m , u>i, one obtains 

u s = sJQi - J s ) , (6.9) 
u m (n-J m ) = 0, (6.10) 
uJt = st/fji-Ji) . (6.11) 

Equation ()6.1()|) gives either u m = or /x = J my leading to two distinct classes of solutions. 
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Figure 6.6: a) Rotational frequency u;(/x) for the four planar bands (marked as A, B, 
C, D) and for the chiral doublet (the vertical line) obtained in the classical model, b) 
Intrinsic-frame trajectory of the rotational frequency vector along those bands. The chiral 
solution corresponds to a straight line perpendicular to the figure plane, and intersecting 
it at the marked point. That perpendicular direction represents u m . 



If uj m = then both oj and J lie in the s-l plane. This gives planar solutions, for 
which the chiral symmetry is not broken. All values of \i are allowed, and the Lagrange 
multiplier must be determined from the length of us calculated in the obvious way from 
()6.9|) and (j6.11|) . Figure l63k shows u versus /i for sample model parameters, extracted 
from the HF PAC solutions with the SLy4 force with no time-odd fields, and listed in 
Table E21 The solutions marked as A and D exist for all values of u, while above some 
threshold frequency, tOthri two more solutions appear, B and C. This threshold frequency 
can be determined by finding the minimum of uo in function of //, and reads 



Wthr 



\Jl — Js\ 



(6.12) 



The value of u t f ir coming from the present HF calculations is rather high, higher than 
lMeV/ft. Since bands B and C are situated far above the yrast line (see Fig. 16 .7|) they 
will not be subject of further analysis. 

For ii=Jmi ah values of io m are allowed, while components in the s-l plane are fixed 
at u s = s s /(J m — J s ) and u>i = si/(J m — Ji). Consequently, the angular momentum has 
non-zero components along all three axes, and the chiral symmetry is broken. For each 
value of u, there are two cases differing by the sign of ou m , and thus giving the chiral 
doublet. The fact that u s and uji are constant leads to the principal conclusion that chiral 
solutions cannot exist for u smaller than the critical frequency 



+ 



Si 



1/2 



(6.13) 



At that frequency, and with u m =0, the chiral solution coincides with the planar band A. 
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Figure 6.7: a) Angular momenta J(u) and b) Energies E(I) for the four planar bands 
and the chiral doublet obtained in the classical model. 



In the u) space, the four planar solutions form a hyperbola in the s-l plane, while the 
chiral doublet corresponds to a straight line perpendicular to that plane. These curves 
are shown in Fig. 16.6b . Figure l6~7k gives the angular momentum in function of rotational 
frequency for all the presented bands. With increasing u, the second moment of inertia, 
dJ/du, asymptotically approaches Ji for bands A and B, and J s for bands C and D. For 
the chiral band, J is exactly proportional to uo with the coefficient J m . Thus, the critical 
spin, J cr it, corresponding to the critical frequency (j6.13J) reads 

J crit \Jm!^crit • (6.14) 

Figure 16.7b summarizes the energies in function of spin. At low angular momenta, the 
yrast line coincides with the planar band D. Then it continues along the planar solution 
A. Since the moment of inertia J m is the largest, beyond the critical frequency the chiral 
solution becomes yrast, thereby yielding good prospects for experimental observation. 

Since the works of Volterra of 1898 |Vol98aj . similar classical particle-rotor models 
have been considered in several contexts different than the chiral rotation. They do 
not always assume a rigid coupling of gyroscopes to the core, but introduce a potential 
interaction between them. Bohr and Mottelson [Boh80a , and later also Kamchatnov 
|Kam90a] examined one gyroscope interacting with an axially-symmetric core, to study 
high-K bands. As far as the magnetic rotation is concerned, Macchiavelli et al. [Mac99aJ 
considered two gyroscopes coupled to a spherical core to investigate the competition 
between the collective rotation and shears mechanism; see Section 12.11 That model was 
generalized by Pasternak |Pas04aj to the case of axially deformed core. 



6.5 Classical estimates of the critical frequency 

A useful feature of the classical model presented in the previous Section is that all its 
parameters can be easily extracted from the standard PAC calculations (see Section ? 
that can be performed by many existing mean-field codes. Then, Eq. (J6.13)) provides 
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an estimate of the critical frequency, and gives insight into where its value comes from. 
Below, the critical frequency and spin are estimated on the basis of the HF PAC results 
of Section 16.31 an d also from the PAC calculations with pairing, performed within the 
Total Routhian Surface (TRS) method. 

The model parameters, J S) J m) Ji, and s s , S[, were obtained from a least-square fit 
of straight lines, Jj = JiUi + s^, i — s,m,l, to the calculated points shown in Fig. 16.41 
Parameter s m was set to zero. Whenever there was a bend in the calculated Ji(u>i) 
dependence (see Section IQj) . the line was fitted in the rotational frequency range below 
the bend. The fitted moments of inertia and initial alignments, together with the resulting 
values of ujf r °£ and I^" t s ( 4 ), are listed in Table IS~2l 

As expected from the irrotational-flow model |Boh75aj . the moment of inertia with 
respect to the medium axis is the largest. Note, however, that the J s : J m : Ji ratio 
is far from that predicted by the irrotational-flow formula (|2.9J) used for the calculated 
triaxiality 7. The inclusion of the time-odd fields G (see Section 13 .4j) increases all the 
moments of inertia, but particularly J~ m , with respect to the case, N, with no time-odd 
fields. This causes a decrease in but the corresponding if^ does not change much, 
because J m is larger. Switching on the fields T results in values of J m between those 
obtained for the cases N and G. The critical frequency always becomes higher than for 
the G fields, and the resulting critical spin is always the highest among all the examined 
sets of time-odd fields. These variations of If°£ are of the order of a few spin units. 

For the chiral bands in 132 La, the full HF TAC solutions were found (see Section IFTTHj) . 
and the obtained values of the critical frequency and spin, u^ r ( t and 1^, are also listed 
in Table IB~2l They are slightly higher than the classical estimates, and 1^, but the 
agreement can be considered good, given the simplicity of the classical model. When the 
time-odd fields are included, values of uj^ r f t and I^ r [ t vary similarly to those of ujf^l and 

Tolas 
crit ' 

To examine the possible role of pairing correlations, PAC calculations within the TRS 
approach were performed for 132 La; see Appendix O for details. Deformation of (3 « 0.20 
and 7 ~ 25° was obtained from energy minimization. In the cranking calculations, TRS 
effectively interchanges values of J s and J7i, and enlarges J m roughly twice with respect to 
the HF results. On the basis of the irrotational-flow model, these changes in the moments 
of inertia can be partly understood as due to the change in triaxiality; the value of 7 
obtained from TRS is closer to 30° and on the other side of 30° than the HF results. 
As governed by Eq. (|6.13jl . the increase in J m significantly lowers the values of uifrit 
and I ciu ^ see Table 16.21 However, inclusion of pairing in the full HF TAC calculations 
(Hartree-Fock-Bogolyubov method) is necessary to verify this trend. 

Values of the critical spin obtained from the HF calculations are rather high as com- 
pared to the spin range in which the supposed chiral bands are observed in nuclei under 
study. A detailed comparison with experiment is given in Section IfHfl 



In the present work, the symbol J denotes the average angular-momentum vector calculated from a 
mean-field solution or the classical angular-momentum vector. The symbol / is used for the spin quantum 
number, as it is measured in experiment. If / is approximated from mean-field or classical results, the 
relation /(/ + 1) = J 2 is used. 
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Table 6.2: The (3 and 7 deformations, parameters of the classical model, J s , m ,i [h 2 /MeV], 
s s j [h], classical estimates for the critical frequency and spin, uf^l [MeV/ft] and lf r ^ [h], 
and full HF TAC results for those quantities, to^ r f t [MeV/h] and I%£ [h], for the iV = 75 
isotopes. HF results with the SLy4 and SkM* forces are listed for the N, G, T time-odd 
fields included. For 132 La, results from TRS calculations are given, too. 
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6.6 Planar HF bands 

The next step in the HF analysis of the vr/i^^ vh~y 2 configuration in the N = 75 isotones 
was to obtain planar solutions, much similar to the shears solutions in 142 Gd; see Sec- 
tion They serve us as starting points in the search for the chiral bands, as described 
in the next two Sections. 

In the N = 75 isotones, planar solutions were found by applying a tilted cranking 
vector to converged non-rotating states. Since in the triaxial minima in 134 Pr and 136 Pm, 

it 



the neutron hu/ 2 states interleave with other negative-parity levels (see Section 
was very difficult to follow the proper configuration diabatically, and the planar bands 
in those minima were obtained only for very low angular frequencies. Figures ID.12MTJT4I 
give sample s.p. Routhians for the planar bands in 132 La. 
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Figure 6.8: Angular momentum vectors of the valence ^11/2 proton particle and neutron 
hole along the HF planar bands found in the N = 75 isotones. Results with no time-odd 
fields are shown; the triaxial solutions in 134 Pr and 136 Pm were calculated with the SkM* 
force, while the remaining results displayed were obtained with SLy4. The angle between 
the two spin vectors is given in each panel. Rotational frequency is given in MeV/h. The 
scale of the plots and the values of u are identical to those in Fig. 15. 3| for the shears band 
in 142 Gd. 

Figure EEB1 shows the angular momenta of the valence hn/2 proton particle and neutron 
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hole for the planar bands in the studied N = 75 isotones. Results with no time-odd 
fields are shown, as an example, because they were found in the largest frequency range, 
but inclusion of the G or T fields very little alters the alignments, like it was already 
demonstrated for the shears band in 142 Gd; see Fig. 15.31 One can see that, indeed, the 
individual spins of the valence nucleons are rather tightly aligned along the short and long 
axis, as it was expected from the PAC calculations in Section 16.31 Interestingly, there is 
almost no difference between the triaxial and oblate minima in 134 Pr and 136 Pm as far as 
this point is concerned. Neither do the four N = 75 isotones differ much in this regard. 
Note, however, that the alignments in the N = 75 isotones are visibly stiffer than for the 
shears band in 142 Gd (see, e.g., the angles between the blades given in Figs. lfi.8l and l5T3~j) . 
This is because the deformation in the considered N = 75 is larger than that in 142 Gd. 

Contrary to the shears band in 142 Gd (see Chapter E"), the deformation is almost 
constant along the planar bands in the N = 75 isotones. The absolute deformation (3 
drops by not more than 0.02 in the frequency range, where the planar solutions were 
found, and the triaxiality varies by about 3°, at most. This is presumably connected to 
the fact that there are only two valence nucleons, compared to four in 142 Gd, and they 
are unable to polarize the core so much. 

The HF planar solutions discussed in this Section start at zero rotational frequency 
(convergence problems with the G and T fields withstanding), and the individual angular 
momenta of the valence particle and hole have positive projections on the u) vector. Thus, 
these solutions have all the characteristics of the planar band A of the classical model; 
see Section 16.41 Indeed, the trajectory of u> along the HF bands almost exactly follows 
the classical hyperbola, as illustrated in Fig. 16.91 Also the calculated energy as function 
of spin agree very well with the classical prediction, see Figs. I6.14H6.1~B1 

6.7 Search for the chiral HF solutions 

Planar HF solutions were easily obtained by applying a small cranking frequency to the 
non- rotating state (at least with no time-odd fields). For chiral bands, analogous task was 
more difficult, because these bands start at a finite frequency, which in the present case is 
not lower than « 0.5 MeV/h. Several level crossings may occur between u — and such 
a high frequency, and it is a priori difficult to identify the required s.p. configuration. To 
follow a given configuration diabatically, it is desirable to find a continuous path linking 
the chiral solution to the non-rotating state. A hint on how to proceed comes from the 
classical model prediction that the chiral band sprouts from the planar solution (at the 
point corresponding to the critical frequency). One can thus restart iterations from the 
converged planar band, rather than from the non-rotating state. Such a way of proceeding 
was adopted in this work, to find the solutions described in the next Section. However, the 
author also found chiral solutions by performing a perturbative search along the planar 
band. That method gives some understanding of why chiral solutions do not appear in 
some cases, and therefore merits a note in this Section. 

Specifically, to each converged point of the planar band, a small additional component, 
u m , of the angular frequency along the medium axis was added. The resulting Routhian 
was diagonalized only once, and the same parity s.p. configuration was required as for the 
planar state. Then, it was checked whether in the resulting (non- self consistent) state the 
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Figure 6.9: Intrinsic-frame trajectory of the angular frequency vector along the HF planar 
bands in the N = 75 isotones, compared to the classical predictions. HF results with the 
SLy4 and SkM* forces are shown for the N, G, T time-odd fields included. 



76 



CHAPTER 6. HARTREE-FOCK CHIRAL RESULTS 



> 

o 



50 



40 



30- i 



— Classical J s /oj s = Ji/u>i 


O HF J s /u s = 


Ji/uji 


— Classical J™ / u) m 


+ HF J™ / CUm 

1 lib I III 




130 Cs 132 La 

i j i i i i i i i i | i 1 i | i 1 i 


134p r 
'I'll 1 1 


136p m 

1 i 1 i 1 i 



20 



3 40 
30~ 



20 



6 



© 
© 



J I I i I I L 



© 



6 
© 



H 1 h 



© 



© 



© 



: | -H+H4f ' - 



+ 



oblate 



<S> 



I 1 I 1 I 1 I 



P+, 



triaxial 



J i I i I i L 



oblate 



I | I | I 



triaxial 



u \ 

:++++ 



J i I i L 



CO 



0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 1 

Rotational frequency [MeV/ft] 



Figure 6.10: The HF and classical-model values of the J m /uj m and J s /uj s w Ji/uji ratios for 
the perturbative search for the HF chiral solutions along the planar bands in the N = 75 
isotones (see text). HF results for the SLy4 and SkM* forces with no time-odd fields are 
shown. 

angular momentum and rotational frequency vectors were parallel or not. The Kerman- 
Onishi theorem [Ker81a requires that their parallelism is a necessary condition for self- 
consistency; see Section EP1 Presuming that in nuclei that are stiff against deformation 
changes, the direction of J is the only degree of freedom, the Kerman-Onishi condition 
is also sufficient. In other words, if J is parallel to u> in the considered state after one 
diagonalization, then it is very probable that further iterations may lead to a converged 
chiral solution. Indeed, in practice it was always the case, and never a chiral solution was 
obtained, in spite of several attempts, if that simple test gave negative result. 
The condition for J and u) being parallel can be written in the form 

Js J-m Jl 

U s UJ m Ui 

Note that the J s /uj s and J\ju)i ratios must be very close to each other in the non- 
selfconsistent state, because the Kerman-Onishi condition is fulfilled for the self-consistent 
planar state. Therefore, the test consists in checking for each point of the planar band if 
Jm/^m is equal to J s /uj s ps Ji/oJi- In fact, this is only reliable if the time-odd fields are 
switched off, because in their presence, the J m /^m ratio calculated perturbatively is sig- 



(6.15) 
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nificantly smaller than the self-consistent result would be 5 . The test was made with u m = 
0.05 MeV/h, and it has been verified that changing this value up to u m rs 0.10 MeV/h 
does not influence the results. The discussed ratios, calculated for all the HF planar bands 
found in the N = 75 isotones, are plotted in Fig. I6.10[ together with the classical-model 
predictions. 

It can be a priori expected that chiral solutions do not appear in the oblate minima 
in 134 Pr and 136 Pm, because of insufficient triaxiality. Indeed, the calculated values of 
the J m /u) m and J s /u> s « Ji/oJi ratios exhibit a complicated behavior, and do not become 
equal one to another at any point. In 130 Cs, as well as in the triaxial minimum in 136 Pm, 
the two ratios clearly approach each other. It seems that the only reason why they do not 
attain equality is that the planar bands were not found up to sufficiently high frequencies, 
because of level crossings. Note, however, that the moment of inertia associated with 
the medium axis, J m /u m , significantly drops with angular frequency, which takes J m /u m 
away from J s /uo s rs Ji/oJu an d defers their equalization to higher frequencies. This does 
not occur in 132 La, where the ratios do become equal, and this happens near the point 
expected from the classical model. As described in the next Section, self-consistent chiral 
solutions were indeed found in this case. 



6.8 Self-consistent chiral Skyrme-HF solutions 

Apart from the perturbative test described in the previous Section, self-consistent chiral 
solutions were sought in the following way. To the mean field of each converged point 
of the planar band, a cranking term was added with the rotational frequency of the 
same magnitude as for the planar solution, but with equal components on all the three 
intrinsic axes. Of course, other choices for u> are possible; the important point is to 
break the chiral symmetry at the starting point by applying a non-zero u m , to render the 
calculations symmetry-unrestricted; see Section IH.HI From the Routhian constructed in 
such a way, regular HF iterations were restarted with the same parity s.p. configuration 
as for the planar state. 

In 130 Cs, 134 Pr, and 136 Pm, as well as for low rotational frequencies in 132 La, the 
iterations converged to planar solutions that differed from the planar starting points 
uniquely by orientation in space (due to different direction of the applied w). This confirms 
in the frame of the Skyrme-HF method that, for the concerned configuration, the chiral 
symmetry is not dynamically broken at low frequencies, as predicted by the classical 
model. In 132 La, for frequencies high enough, converged solutions were obtained with the 
angular momentum having non-zero components on all the three intrinsic axes, and thus 
violating the chiral symmetry, R T . These are the first fully self-consistent results that 
confirm the spontaneous breaking of this symmetry in rotating atomic nuclei. 

The found fragments of chiral bands were then used as starting points to obtain chiral 
solutions for lower and higher angular frequencies. Calculations were performed with a 
frequency step of 0.02 MeV/fo. At a certain value of decreasing u, the planar orientation 
of the angular momentum was regained in the intrinsic frame, and the solution merged 

5 This can be seen, e.g., from the comparison of non-selfconsistent and self-consistent PAC results. The 
reason is that the relevant components of the time-odd fields only become active after self-consistency is 
achieved for a non-zero cranking frequency. 
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into the previously obtained planar band 6 . This confirms the classical prediction that 
the planar and chiral bands have a common point. In a natural way, that junction value 
of us can be regarded as the Skyrme-HF result for the critical frequency, and is denoted 
in the following as uj^,^. From Figs. ITTT21 to ITTul it can be seen that the chiral and 
planar s.p. Routhians do indeed coincide at u — oj^ r f t . Values of to^-ft are collected in 
Table 16.21 and are discussed in the next Section. For highest frequencies, chiral solutions 
were obtained up to a certain value of u, and all attempts to go higher caused iterations to 
fall into another, although also chiral, minimum. This is probably due to multiple smooth 
crossings of occupied and empty levels, particularly in neutrons; see the s.p. Routhians in 
Figs. EH to EH That other minimum is not build on the vr/i^^ vh x y 2 configuration, 
and will not be subject of further studies here. Although the chiral solutions have been 
found in a rather narrow u interval, of about 0.1 MeV /h, the accompanying increase in 
uj m is significant, from zero to about 0.4 MeV/h; see Fig. 16. Ill 
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Figure 6.11: Trajectory of the angular frequency vector along the HF chiral bands in 
132 La, compared to the classical predictions. Projections onto the s-m and l-m planes of 
the intrinsic frame are shown. The Figure gives the HF results with the SLy4 and SkM* 
forces for the N, G, T time-odd fields included. 

Similarly to the planar solutions described in the previous Section, there is almost 

6 With the exception of the SLy4 case with the T time-odd fields. Although a chiral solution has been 
found for this force, it could not be connected to the planar band, and it is not sure whether it is built 
on the 77^^/2 v ^ l ~\\i2 configuration. In the following, results for this band are not discussed. 
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Figure 6.12: Alignments of the valence /iii/2 proton particle on the short, medium and 
long intrinsic axis along the HF chiral band in 132 La. HF results with the SLy4 and SkM* 
forces are shown for the N, G, T time-odd fields included. 



no deformation change along the HF chiral bands in 132 La. This yields good prospects 
for the validity of the classical descriptions of those bands; see Section 16.41 Figure 16.111 
gives the intrinsic-frame trajectories of the angular frequency vector along the HF chiral 
solutions. In each case, the HF trajectory is almost a straight line, parallel to the medium 
axis, as predicted by the classical model. The only difference is that u)^ r ( t is a bit higher 
than cj^it, and the line is shifted along the planar band to higher frequencies. 

In the vicinity of uj^[ t , the hwn Routhian associated with the valence neutron hole 
enters a region of high density of negative-parity levels interacting among themselves. 
Therefore, it is doubtful whether the considered hole can still be identified with a single 
s.p. state. However, the hn/ 2 proton particle is still well separated. Figure 16.121 gives 
its spin alignments on the three intrinsic axes along the chiral band. The plot confirms 
the stiff character of those alignments. Like in the case of u), only the component on the 
medium axis changes, the other two remaining nearly constant. 



6.9 Comparison with experiment 

Unfortunately, mean-field results for the chiral rotation are not directly comparable with 
experimental data, as those of the PRM. This is because the pure mean-field approach 
does not take into account the interaction between the left- and right-handed solutions, 
which are exactly degenerate; see Section l2~2l Thus, splitting between the chiral partners 
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[E(I + 1) — E(I — l)]/2 for the candidate chiral bands in the four N = 75 isotones. Solid 
and dashed lines connect the points corresponding to odd and even spins, respectively. 



cannot be reproduced, nor their exact energies. It is argued in the literature that the 
mean-field chiral solution can be viewed as a kind of average of the two partner bands, 
and one can thus compare the mean trends. One can also speculate about the value of the 
critical frequency. It should correspond to the kink in the I (no) dependence, as predicted 
by the PRM and by the classical model; see Figs. 12.51 and 16.71 respectively. However, 
there seems to be no distinct kinks in the experimental bands in the N = 75 isotones; see 
Fig. EI31 

Figures I6.14H6.161 give a comparison of calculated and experimental energies. At low 
spins, where the supposed chiral partners have not been observed, the level energies are 
reasonably reproduced by the HF planar solutions. It may mean, on the one hand, that 
the yrast band is not yet perturbed by interactions with the chiral partner, and, on the 
other hand, that the HF moments of inertia are correct. Roughly at the spin where the 
chiral partners commence to be visible, the HF planar solutions cease following the yrast 
bands, which significantly change their behavior. One may argue, therefore, that this is 
due to the interaction between the left- and right-handed minima. 

In 132 La, the HF results for the critical spin, I^ r [ t = 15.9 — 18.5h (see Table IB~2l or 
Fig. I6.15|) are rather high as compared to the spin range, in which the two candidate chiral 
partners are observed. Yet, the classical estimate, If^ = 9.2h, evaluated for the TRS 
PAC results with pairing (see Section 16. 5|) is already below that range. It means that the 
inclusion of pairing in the calculations may be important for the correct interpretation of 
the data. But it also suggests that the actual critical spin is somewhere between the HF 
and TRS results, that is roughly in the middle of the spin range of the chiral partners. 
The HF-PAC-based values of lf r a ^ calculated for 130 Cs and 136 Pm also lie at half-spin- 
range of the partner bands; see Figs. 16.141 and 16. lTfl The closeness of the side bands to the 
critical spin may imply that those structures represent the transition from planar to chiral 
rotation, as discussed in Section l2~2l Although this transition is abrupt in the cranking 
model, and even in the PRM as indicated by the sharp kink in the I(uj) dependence in 
Fig. 12.5b . it may be rather smooth and complex in real nuclei. This is an interesting topic 
for study including techniques beyond the mean field. 
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Figure 6.14: Energies from the HF TAC calculations and from the classical model, com- 
pared with the experimental data on the candidate chiral partners in 130 Cs. HF results 
with the SLy4 and SkM* forces are shown for the N, G, T time-odd fields included. 
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Figure 6.16: Same as Fig. I6.14| but for the bands in 134 Pr and 136 Pm. 
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Chapter 7 

Conclusions and outlook 



To conclude, a numerical tool has been constructed that allows for Skyrme-Hartree- 
Fock (Skyrme-HF) Tilted- Axis Cranking (TAC) calculations with no symmetries imposed. 
Thus, it can be used for investigating the spontaneous breaking of various symmetries in 
rotating nuclei. As a first application, the phenomena of the magnetic and chiral rotation 
were studied. For their cranking description, both two-dimensional (planar) and three- 
dimensional (chiral) TAC solutions were obtained, which is a robust test of the method. 
The present results provide one of the first fully self-consistent proofs for the possibility of 
nuclear rotation about an axis tilted with respect to the principal axes. They constitute 
the first fully self-consistent description of the chiral rotation, and were preceded only by 
a single Relativistic Mean Field calculation [MadOOa] . as far as the magnetic rotation is 
concerned. 

Before going to the full HF TAC calculations, properties of the valence nucleons 
were examined within the PAC frame and on the basis of pure symmetry considerations. 
It was found that the valence particles and holes align their angular momenta along 
the short and long axis of the nucleus, respectively. It was argued that their spins are 
bound to those axes rather tightly, and possibly more tightly in triaxial nuclei than in 
axially-symmetric ones. By considering a very simple model it was concluded that pairing 
correlations may soften those alignments. 

For the study of the shears bands, the nh^^ uh^, 2 configuration in 142 Gd was chosen. 
A substantially triaxial solution was obtained. The HF results corroborate that an impor- 
tant portion of the angular momentum is generated by the shears mechanism, as expected 
from earlier investigations. It was found that the aligning valence nucleons polarize the 
core quite strongly, and the deformation decreases with spin. However, the rotating core 
also generates angular momentum, and actually wins the competition with the shears 
mechanism. This may be due to overestimated deformation or the lack of pairing in the 
calculations. It is also possible that triaxiality acts against the shears mechanism by stiff- 
ening the alignments of the valence nucleons; such an effect has not been considered so 
far. The too big collective moment of inertia causes that agreement with experiment is 
not satisfactory. 

The chiral solutions were sought in four N = 75 isotones, 130 Cs, 132 La, 134 Pr, 136 Pm, 
with the fth\y 2 vh~y 2 configuration. Strongly triaxial solutions were found, but nearly 
oblate minima proved to be lower in energy for 134 Pr and 136 Pm. One-dimensional crank- 
ing about the three axes premised that the alignments of the valence nucleons are very 
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stiff and that the core responds to rotation linearly, like a rigid body. On the basis of 
these results, a classical model was constructed, in which the chiral rotor is represented 
by two gyroscopes coupled to a triaxial rigid body. The solutions to that model agree 
surprisingly well with the full HF TAC results, and elucidate the salient features of the 
chiral rotation. The HF chiral solutions were found only in 132 La. It was established that 
the chiral rotation appears only above a certain critical value of the angular frequency; a 
conclusion that had been known before as a purely numerical result. In the present work, 
the origin of the critical frequency is explained in an illustrative way within the classical 
model, that also gives an analytical estimate for that quantity. The HF chiral solutions 
correctly reproduce the mean E(I) behavior of the observed doublet bands in 132 La. The 
calculated critical frequency is rather high as compared to the spin range in which the 
chiral partner bands are observed. PAC calculations with pairing suggested that pairing 
correlations may reduce its value. It is also possible that the bands in the N = 75 isotones 
represent the transition from planar to chiral rotation. 

One purpose of the present study was to investigate the influence of the Skyrme time- 
odd fields on the magnetic and chiral results. The inclusion of those fields changes mainly 
the moments of inertia, as expected, and causes convergence problems, mainly at low 
cranking frequencies. It does not alter, however, the qualitative features of the solutions. 

Granted, the pairing correlations do influence the shears-rotor competition in the 
magnetic bands and the critical frequency in the chiral rotation. Their inclusion in the 
HF TAC calculations in the frame of the Hartree-Fock-Bogolyubov method would be very 
desirable. The pure mean-field approach is not capable of estimating the splitting between 
the chiral partner bands, because it does not take into account the interaction between 
the left- and right-handed solutions. Inclusion of that interaction requires techniques 
beyond the mean field, preferably the restoration of broken chiral symmetry. To describe 
the smooth transition from planar to chiral rotation, it should be the projection before 
variation. 



Appendix A 

Invariance of the time-even density 



This Appendix gives the proof that unitary mixing of states in a Kramers pair does not 
alter the time-even part of the density matrix, even if only one particle is put on such a 
pair. This lemma was used in Section [Ql 

The density matrix, p, of an arbitrary Slater determinant, is defined as 

p af3 = (#|a+a a |^) . (A.l) 

It can always be expressed as a sum of contributions from the single-particle (s.p.) states, 
\vi), of which |\&) is constructed, 

Pap = ^(0|a^a+a Q a+|0) . (A. 2) 

i 

The density p can be uniquely decomposed into a time-even part, p( + \ and a time-odd 
part, p(~\ 

p = p (+) + p(") , f V+)f = +p« , T + p ( -)T = -p(-) , (A.3) 
in the following way: 

P (+) = ^(p + P), P { - ] = \(P-P)i (A.4) 

where p is the time-reversed density, i.e. the density of the state T|\I>), 

p = f + pf , p af3 = ((V\f + )a+a Q (f I*)) . (A.5) 

Since the time-reversal operator, T, is anti-unitary and anti-hermitian when acting on s.p. 
states, one can always choose the s.p. basis so that it consist of pairs of states, (|p), |/2)), 
mutually reversed in time up to phase factors s M , (Kramers pairs), 

T\p) = s^lp.) , f\p) = Sp ,\p) , = -s M . (A.6) 

In such a basis, the time-reversed density matrix, p, is related to the original density, p, 
by a simple formula, 

Pa/3 = S a S*pP* af3 . (A.7) 
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See, e.g., jDob87a] for detailed derivations. 

If the s.p. Hamiltonian, h, is even under transformation through the time-reversal, 

f + hf = +h , (A.8) 

then its eigenstates can be chosen so that they form a basis like above. The states \p) 
and \p) have equal energies (Kramers degeneracy). Suppose now that of such eigenstates 
one constructs a Slater determinant, and in at least one Kramers pair only one state 
is occupied, say \p). This may happen if the number of particles is odd or if there are 
particle- hole excitations. According to (jA.2|) . the contribution to the density matrix, p, 
from the particle on \p) is 

Ap Q/3 = (0|a M aJa a aJ|0) = 5^5^ , (A. 9) 

where there is no summation over p, and p has been expressed in the basis of the same 
s.p. states, of which is built. By using (|A.4|) and (|A.7j) one finds that such a single 
particle contributes both to the time-even and to the time-odd density, 

A p ( af3 = ^(<W<W + <WW) > A P ( a(3 = ^(<W<W - <W<W) • ( A ' 10 ) 

Since the states \p), \p) are degenerate, all their unitary combinations are equally good 
eigenstates of h and can serve for constructing Slater determinants. The theorem used in 
Section IP1 states that if the particle is put not on \p), but on a state |//}, 

\li')=p\li) + q\P) , p*p + q*q=l, (A.ll) 

then its contribution to the time-even density does not change. Indeed, its new contribu- 
tion, Ap', to the total density equals 

Va/3 = (0\(a^p* + a Ii q*)a+a a (pa+ + qa+)\0) (A.12) 
= P>W<W + P*q S afi8/3 fl + q*p5 af ,5p [1 + q*qS a p5p li , (A. 13) 

where the same basis has been used as in ()A.9|) . Again, from (jA.4|) and (jA.7|) one calculates 
directly that Ap'^ = Ap( + \ and the time-odd term, Ap(~\ changes. Note yet, from 
flA.lOjl . that if both states are occupied in a Kramers pair, then their contributions to the 
time-odd density cancel. 



Appendix B 

Alignment and decoupling vectors 



For a .Dj-symmetric single-particle (s.p.) Hamiltonian, the alignment and decoupling vec- 
tors of eigenstates forming a Kramers pair are expressed in terms of three real parameters. 
Thus, the Equations (|4.H| I4.4| 14. 5j) of Section 14. HI are derived. 

Consider a Kramers pair, |/2)), as defined by (jA.6|) . For the state one can 

define the real alignment vector, «J M , and the complex decoupling vector, D^, 

J"={ti\y\n), D» = W\t). (B.l) 

Although the decoupling vector changes its phase when \p) change their phases, the 
relative phases of its components do not depend on the phase convention. Since the 
angular-momentum operator is odd under the time reversal, it can easily be verified that 

J" = (fl\ J\fi) = -J» , D* = (fi\J\ri = D»* (B.2) 

Here, we are concerned by properties of a single Kramers pair formed by eigenstates 
of a s.p. Hamiltonian, h, that is even under the D\ group. Therefore, it is possible to 
choose the states of the pair as eigenstates of either of the three signatures, Ri, where 
i = x,y, z, but only one at a time, because the signature operators do not commute 
among themselves. This results, respectively, in three formally different pairs, (\fii), 
that correspond to just three different bases in the same two-dimensional eigenspace of 
h. Here, it is also convened that the states have eigenvalues of — i under the action 
of Ri, while the eigenvalues of are +i. 

The fact that and are eigenstates of Ri, together with the transformation rules 
of the components, Jj, of the angular momentum operator under the three signatures, 

RtJA = { !f j -Z ) > (B.3) 

induces limitations on the components, and D^, of the alignment and decoupling 
vectors. Namely, 

Jil% = | non-zero j = i ^ = ( j = i 

i 1 j ^ i ' J 1 non-zero j ^ i 

In other words, is confined to the axis i and -D M,; to the plane perpendicular to that 
axis. 
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In order to establish connections between the alignment and decoupling vectors as- 
sociated with the three pairs, one needs to express the states of one pair as 
linear combinations of those in the remaining two. Below, states \n x ), \p, x ) and \/J. y ), \p. y ) 
are expressed as linear combinations of \pb z ), \fi z ) by diagonalization of R x and R y in the 
basis of \fi z ), \fi z ), after having found the corresponding two-by-two matrices. Taking into 
account that R x and R y are odd under R z , 



R z R X R Z 



-R x 



R z R y Rz Ry i 



(B.5) 



one arrives at a conclusion that the diagonal matrix elements of R x and R y in the \z), \z) 
basis vanish. The off-diagonal ones depend on the phases of \z), \z), but are restricted 
by the fact that the signatures are unitary anti-hermitian operators. Moreover, the off- 
diagonal elements of R x and R y are connected by the equality 



(z\R x \z) = i(z\Ry\z) , 
which is due to the multiplication rule 

Rx R y Rz • 



(B.6) 



(B.7) 



All this means that, with a proper choice of phases in \{i z ), |/i 2 ), one can always make 
that the corresponding matrices of Ri be related to the Pauli matrices 1 , <7j, 



Ri = -ibi 



(B.8) 



By diagonalization of R x and R y one obtains, with an arbitrary choice of phases: 



IfO = y y(|/^) + 1/4*)) , 



If 1 *) = - 



(B.9) 



~(\Hz) + i\^z)) , \Hy) 



—(i\li x ) + |//*)) • 



;b.io) 



These formulae, together with relations (|B.2)1 . allow to express J^ x , D^ x and J^ v , D 1 *" in 
terms of J^ z , . Taking into account the constraining conditions (|B.4|) on the left-hand 
side, one obtains the following set of equations: 



( J«° , 0, 0) = (Re Re 0) , (0, Df, D?) = (-Im D 



X ) 



-Im Dp, 



;b.ii) 



(0, J^, 0) = (-Im Dp, -Im Dp, 0) , (£>£», 0, = (-zRe -*Re D^, ) . 
(B.12) 

lr rhis assures compatibility of notation with the case of a particle possessing only the intrinsic spin and 
no space degrees of freedom, where the general definition, Ri — exp(— in Jj), leads explicite to Ri = —i&i. 
For such a particle, the states of spin | f) and | .[) have the z-signature signature eigenvalues of — i and 
+i, respectively, and play the role of \/j, z ) and \£L Z )- 
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A glance at them convinces us that all of the quantities Dj' express through the 
three "diagonal" components, Jf\ One eventually arrives at: 

J"" = ( Jg*, 0, 0) , D»* = (0, J^, -U?) , (B.13) 

J»y = (0, J*<, 0) , = (-jjf,0, Jf ) , (B.14) 

J"« = (0,0, Jf) , D»* = (J^,-U^,0) . (B.15) 

If the Hamiltonian, h, has no extra symmetries, then can assume arbitrary values. 
If h is axially symmetric, say with respect to the z axis, then the states |// z ), \p, z ) are 
eigenstates of J z , which leads to quantization of J^ z . In fact, J^ z = +1/2,-3/2,..., 
because R z = exp(—inJ z ), while in the adopted convention R z \fi z ) = —i\fjL z ). For states 
\fi x ) and [fly), defined by ()B.9j) and (jB.lOjl . one easily finds 

J? = ^Re(fx z \ J + + L\fi z ) , = ~Re</^| J+ - J-\ft z ) , (B.16) 

where J+ = J x + iJ y and J_ = J x — i.J y are the ladder operators, that increment and 
decrement the magnetic quantum number, J z , of an eigenstate, | J z ), of J z , 

J+\J Z )~\J Z + 1) , J_\J Z )~\J Z -1) . (B.17) 

One can see, therefore, that the matrix elements in (|B.16|) can be non-zero only if |/i z ) 
and \p z ) differ in J z by one, that is if J^ z = 1/2. In such a case, (jj, z \J_\p, z ) = 0, and 
J%* = J*f*. These results can be summarized as 

, f (Jf\Jf, J?) for J,; 11 = 1/2 , 
W.^^) = {; ,0,V;) " for jj" = 3/2, 5/2, ... < R18 > 



The parameter J^ x is not restricted by the above kinematic conditions, and is usually 
called decoupling parameter. 
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Appendix C 

Technical details of the calculations 



This Appendix summarizes all technical details of the present Hartree-Fock (HF) calcu- 
lations, and of the auxiliary calculations performed in Chapters El and |B1 within the Total 
Routhian Surface (TRS) approach. 

All the HF calculations presented in this work were performed with the help of the 
code hfodd, version (v2.05c), which is described in Section Twelve spherical shells of 
the Harmonic Oscillator were taken as the single-particle (s.p.) basis. It has been verified 
for the one- dimensional cranking in 132 La that varying the number of shells from 10 up to 
16 changes the quantities important for the present study, like deformation, alignments, 
moments of inertia, by less than 1%. The calculations were performed with two Skyrme 
parameter sets, SLy4 [Cha97aJ and SkM* [Bar82aJ. Attempts were undertaken to obtain 
each solution with three sets of the HF time-odd fields |Dob95a[ |Dob97a included, as 
listed in Table 1X2"! and explained in Chapter 13.41 Those sets are denoted in the text as N, 
G, T. Since the treatment of the pairing correlations is not yet implemented in hfodd 
in the Tilted-Axis Cranking (TAC) mode, pairing was not included. In all Principal-Axis 
Cranking (PAC) calculations, the P, R y , Sj, Sj symmetries (and their products) were 
imposed, while only P was kept in the TAC case. The respective parity-signature and 
parity configurations (see Section l3"l)j) for all the considered solutions are given in Tabs. EH) 
and l6.ll The configurations were followed diabatically in each band, either by choosing the 
occupied state by hand, whenever a crossing was about to occur, or by using the diabatic- 
blocking techniques described in Section 13.61 In all planar solutions, the direction of 
u) was chosen so that it had equal components on two axes of the program frame (see 
Section 13 .3J) and zero component on the third axis. In the chiral bands, uj had equal 
projections on all the three program axes. In each case, the direction of uj was kept fixed, 
and only its magnitude was changed. The axial and planar bands were calculated with 
an uj step of 0.05 MeV/fr, while the step of 0.02 MeV/h was used for the chiral solutions. 
With no time-odd fields included, the iterations for each value of u were restarted from 
the previously converged solution at the prior value of uj. Then, each point converged 
with no time-odd fields served as a starting point for calculations with the time-odd fields 
included. As explained in Section l3~Bl the convergence rate is sometimes very low in TAC 
calculations. To deal with this difficulty, all solutions were first obtained at the level 
of convergence defined by the hfodd flag EPSITE=0. 000001; see the hfodd manual in 
|Dob00aj . Afterwards, additional iterations were repeated in portions until the functional 
energy ceased changing from one iteration to another by more than 1 eV. This criterion 
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APPENDIX C. TECHNICAL DETAILS OF THE CALCULATIONS 



could be fulfilled after hundreds of iterations at higher frequencies, uj > 0.4 MeV/h, and 
after thousands of iterations for lower frequencies, u < 0.4 MeV/h. This process could 
have been much faster, had the technique of resetting uj parallel to J been known earlier; 
see Section 13.61 The constraint on u> x J was not used, because it turned out to be very 
ineffective; see Section l3~Hl 

To estimate the possible influence of the pairing correlations, auxiliary calculations 
within the TRS method were performed in the present work. In Chapter El the energy 
surface for the ground state in 142 Gd was examined, and one- dimensional cranking for 
132 La was considered in Chapter El The calculations were done by using the TRS software 
developed by Satula and Wyss |Sat94aj . The method employs the phenomenological 
Woods-Saxon potential and a correction for the liquid-drop energy. Self-consistency is 
assured in the pairing channel, and the strength of the monopole pairing is adjusted to 
experimental data following the average gap method. The quadrupole pairing is also 
included, with the coupling constant chosen so that the pairing field possess the local 
galilean invariance. Approximate projection onto good particle number is implemented 
within the Lipkin-Nogami method. The code does not allow for TAC calculations. See 
|Sat94aj for all details and further references. 



Appendix D 

Single-particle Routhians 



In this Appendix, all Figures showing the Hartree-Fock (HF) single-particle (s.p.) Routhi- 
ans are collected. 

The Routhians marked as Axial in the Figures' captions come from one-dimensional 
cranking, with imposed parity and signature. The Routhians denoted Planar and Chi- 
ral were obtained from Tilted-Axis Cranking (TAC) calculations, with only the parity 
imposed. In both kinds of Routhians, solid lines are used for positive-parity levels, and 
dashed lines for negative-parity levels. In the case of Principal- Axis Cranking (PAC), 
levels of positive signature are plotted in black, and levels of negative signature (—i) 
in gray. Each s.p. level is labeled with its total ordinal number and its ordinal number in 
the given symmetry block. For the PAC or TAC Routhians, the second number refers to 
the parity-signature or parity blocks, respectively. For example, the state marked as 62 
14 — h is the 62nd state among all and the 14th of negative parity and positive signature. 
The last item in each label indicates whether the given state is occupied (1) or not (0). 
This information is not provided for states far from the Fermi energy. For clarity, the 
Routhians bearing the hu/ 2 nucleons are marked with points, put at the values of the 
cranking frequency, at which the HF solutions were obtained. The proton particles and 
the neutron holes are marked with full and open circles, respectively. The symbols N, G, 
T, used in the captions, refer to calculations with various Skyrme time-odd fields included; 
see Section EOl The labels oblate and triaxial, used for 134 Pr and 136 Pm, denote the nearly 
oblate and triaxial minima found in those nuclei, see Section IB~2l Figures TP . 1 2HD . 14l show 
the Routhians for the planar and chiral bands in 132 La in the same plots. A vertical line 
is drawn at the value of the HF critical frequency, ui^ r f t ; see Section I6~H1 The Routhians 
left to that line correspond to the planar band, and the Routhians to the right correspond 
to the chiral band. 
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Short axis 



Long axis 
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Figure D.l: Axial Routhians in 142 Gd for SLy4 force with N time-odd fields. 
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Figure D.2: Planar Routhians in 142 Gd with N time-odd fields. 
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Figure D.3: Planar Routhians in 142 Gd with G time-odd fields. 
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Figure D.4: Planar Routhians in 142 Gd with T time-odd fields. 
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Figure D.5: Axial Routhians in 130 Cs for SLy4 force with N time-odd fields. 
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Figure D.8: Axial Routhians in oblate 134 Pr for SLy4 force with iV time-odd fields. 
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Figure D.9: Axial Routhians in triaxial 134 Pr for SkM* force with iV time-odd fields. 
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Figure D.10: Axial Routhians in oblate 136 Pm for SLy4 force with N time-odd fields. 
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Figure D.ll: Axial Routhians in triaxial 136 Pm for SkM* force with N time-odd fields. 
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Figure D.12: Planar and Chiral Routhians in 132 La with N time-odd fields. 
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Figure D.13: Planar and Chiral Routhians in 132 La with G time-odd fields. 
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Figure D.14: Planar and Chiral Routhians in 132 La with T time-odd fields. 



APPENDIX D. SINGLE-PARTICLE ROUTHIANS 



Bibliography 



[Abo92a] Y. Aboussir et al, Nucl. Phys. A549, 155 (1992) 

[AmiOla] Amita, A.K. Jain, B. Singh, Table of magnetic dipole rotational bands 
(www.nndc.bnl.gov, 2001) 

[Ant02a] M.S. Antony, Nuclide Chart 2002 (Strasbourg, 2002) 

[Bal92a] G. Baldsiefen et al, Phys. Lett. B275, 252 (1992) 

[Bal94a] G. Baldsiefen et al, Nucl. Phys. A574, 521 (1994) 

[Bar82a] J. Bartel et al, Nucl. Phys. A386, 79 (1982) 

[BarOla] R.A. Bark et al, Nucl. Phys. A691, 577 (2001) 

[BeaOla] C.W. Beausang et al, Nucl. Phys. A682, 394c (2001) 

[Ben02a] M. Bender et al, Phys. Rev. C65, 054322 (2002) 

[Boh53a] A. Bohr, B.R. Mottelson, Mat. Fys. Medd. Kong. Dansk. Vid. Selsk. 27, No. 
16 (1953) 

[Boh75a] A. Bohr, B.R. Mottelson, Nuclear Structure (Benjamin, New York, 1975) 
[Boh80a] A. Bohr, B.R. Mottelson, Phys. Scripta 22, 461 (1980) 
[Cha97a] E. Chabanat et al, Nucl. Phys. A627, 710 (1997) 
[Chm97a] S. Chmel et al, Phys. Rev. Lett. 79, 2002 (1997) 
[Cla92a] R.M. Clark et al, Phys. Lett. B275, 247 (1992) 

[ClaOOa] R.M. Clark, A.O. Macchiavelli, Annu. Rev. Nucl. Part. Sci. 50, 1 (2000) 
[Dav58a] A.S. Davydov, B.F. Fillipov, Nucl. Phys. 8, 237 (1958) 
[DimOOa] V.I. Dimitrov, S. Frauendorf, F. Donau, Phys. Rev. Lett. 84, 5732 (2000) 
[DimOOb] V.I. Dimitrov, F. Ddnau, S. Frauendorf, Phys. Rev. C62, 24315 (2000) 
[Dim02a] V.I. Dimitrov, F. Donau, S. Frauendorf, |nucl-th/02li"o"63"V l (2002) 
[Dob84a] J. Dobaczewski, H. Flocard, J. Treiner, Nucl. Phys. A422, 103 (1984) 

111 



112 



BIBLIOGRAPHY 



[Dob87a] J. Dobaczewski, Wybrane zagadnienia teorii jadra atomowego (skrypt do 
wykladu) (1987) 

[Dob95a] J. Dobaczewski, J. Dudek, Phys. Rev. C52, 1827 (1995) 

[Dob97a] J. Dobaczewski, J. Dudek, Phys. Rev. C55, 3177(E) (1997) 

[DobOOa] J. Dobaczewski, J. Dudek, Comp. Phys. Comm. 102, 166 (1997); 102, 183 
(1997); 131, 164 (2000) 

[DobOOb] J. Dobaczewski et al, Phys. Rev. C62, 014310 (2000); C62, 014311 (2000) 

[Dob04a] J. Dobaczewski, P. Olbratowski, Comp. Phys. Comm. 158, 158 (2004) 

[Dud04a] J. Dudek, N. Schunck, private communication (2004) 

[Eng75a] Y.M. Engel et al, Nucle. Phys. A249, 215 (1975) 

[Eng03a] J. Engel et al, Phys. Rev. C68, 025501 (2003) 

[Fra91a] S. Frauendorf, T. Bengtsson, AIP Conf. Proc. 259, 223 (1991) 

[Fra92a] S. Frauendorf, F.R. May, Proc. of the Int. Conf. on Nuclear Structure at High 
Angular Momentum (Ottawa, 1992) 

[Fra93a] S. Frauendorf, Nucl. Phys. A557, 259c (1993) 

[Fra96a] S. Frauendorf, J. Reif, G. Winter, Nucl. Phys. A601, 41 (1996) 

[Fra96b] S. Frauendorf, J. Meng, Zeit. fur Phys. A356, 263 (1996) 

[Fra97a] S. Frauendorf, J. Meng, Nucl. Phys. A617, 131 (1997) 

[Fra97b] S. Frauendorf, J. Reif, Nucl. Phys. A621, 736 (1997) 

[Fra97c] S. Frauendorf, Zeit. fur Phys. A358, 163 (1997) 

[FraOOa] S. Frauendorf, Nucl. Phys. A677, 115 (2000) 

[FraOla] S. Frauendorf, Rev. Mod. Phys. 73, 463 (2001) 

[FraOlb] S. Frauendorf, Acta Phys. Pol. 32, 2661 (2001) 

[Fra03a] S. Frauendorf, private communication (2003) 

[Fri87a] H. Frisk, R. Bengtsson, Phys. Lett. B196, 14 (1987) 

[Gol53a] H. Goldstein, Classical Mechanics (Addison- Wesley, Cambridge, 1953) 

[Gro04a] E. Grodner et al, Int. Journ. of Mod. Phys. E13, 243 (2004) 

[Gro04b] E. Grodner, J. Srebrny, private communication (2004) 

[HarOla] D.J. Hartley et al, Phys. Rev. C64, 31304 (2001) 



BIBLIOGRAPHY 

[HecOla] A.A. Hecht et al, Phys. Rev. C63, 51302 (2001) 

[Kam90a] A.M. Kamchatnov, J. Phys. G16, 1203 (1990) 

[Ker81a] A.K. Kerman, N. Onishi, Nucl. Phys. A361, 179 (1981) 

[KoiOla] T. Koike et al, Phys. Rev. C63, 61304 (2001) 

[Koi03a] T. Koike et al, Phys. Rev. C67, 44319 (2003) 

[Koi03b] T. Koike et al, AIP Conf. Proc. 656, 160 (2003) 

[Kuh92a] A. Kuhnert et al, Phys. Rev. C46, 133 (1992) 

[Lal96a] G.A. Lalazissis, M.M. Sharma, P. Ring, Nucl. Phys. A597, 35 (1996) 

[Li02a] X.F. Li et al, Chin. Phys. Lett. 19, 1779 (2002) 

[Lie02a] R.M. Lieder et al, Eur. Phys. Journ. A13, 297 (2002) 

[Mac98a] A.O. Macchiavelli et al, Phys. Rev. C57, R1073 (1998) 

[Mac98b] A.O. Macchiavelli et al, Phys. Rev. C58, R621 (1998) 

[Mac98c] A.O. Macchiavelli et al, Phys. Rev. C58, 3746 (1998) 

[Mac99a] A.O. Macchiavelli et al, Phys. Lett. B450, 1 (1999) 

[MadOOa] H. Madokoro et al, Phys. Rev. C62, 061301 (2000) 

[Mer02a] E. Mergel et al, Eur. Phys. Journ. A15, 417 (2002) 

[Mol95a] P. Mdller et al, At. Dat. Nucl. Dat. Tabl. 59, 185 (1995) 

[Nil55a] S.G. Nilsson, Mat. Fys. Medd. Dansk. Vid. Selsk. 29, No. 16 (1955) 

[Olb02a] P. Olbratowski et al, Acta Phys. Pol. B33, 389 (2002) 

[Olb02b] P. Olbratowski et al, inucl-th/0211075| (2002) 



[Olb04a] P. Olbratowski et al, accepted in Phys. Rev. Lett. 

[Pas04a] A.A. Pasternak, private communication (2004) 

[Pen03a] J. Peng, J. Meng, S.Q. Zhang, Phys. Rev. C68, 044324 (2003) 

[Pen03b] J. Peng, J. Meng, S.Q. Zhang, Chin. Phys. Lett. 20, 1223 (2003) 

[Pet96a] CM. Petrache et al, Nucl. Phys. A597, 106 (1996) 

[Rai03a] G. Rainovski et al, Phys. Rev. C68, 024318 (2003) 

[Rai03b] G. Rainovski et al, Journ. Phys. G29, 2763 (2003) 



114 BIBLIOGRAPHY 

[Rei95a] P.G. Reinhard, H. Flocard, Nucl. Phys. A584, 467 (1995) 
[Rie60a] B. Riemann, Abh. K6n. Ges. Wiss. 9, 1 (1860) 

[RinOOa] P. Ring, P. Schuck, The Nuclear Many-Body Problem (Springer, Berlin, 2000) 

[Rob03a] S.P Roberts et al, Phys. Rev. C67, 057301 (2003) 

[RzaOla] T. Rzaca-Urban, Acta Phys. Pol. B32, 2645 (2001) 

[Sat94a] W. Satula, R. Wyss, Phys. Rev. C50, 2888 (1994) 

[Sta88a] W. Starzecki et al, Phys. Lett. B200, 419 (1988) 

[StaOla] K. Starosta et al, Nucl. Phys. A682, 375c (2001) 

[StaOlb] K. Starosta et al, Phys. Rev. Lett. 86, 971 (2001) 

[Sta02a] K. Starosta et al, Phys. Rev. C65, 044328 (2002) 

[Sta03a] K. Starosta, P. Joshi, private communication 

[Sug97a] M. Sugawara et al, Zeit. fur Phys. A358, 1 (1997) 

[Tho62a] D.J. Thouless, J.G. Valatin, Nucl. Phys. 31, 211 (1962) 

[Tim03a] J. Timar et al, Eur. Journ. of Phys. A16, 1 (2003) 

[Uma91a] A.S. Umar et al, Phys. Rev. C44, 2512 (1991) 

[Vam04a] C. Vaman et al, Phys. Rev. Lett. 92, 032501 (2004) 

[Vol98a] V. Volterra, Acta Mathematica 22, 201 (1898) 

[Woo54a] R.D. Woods, D.S. Saxon, Phys. Rev. 95, 577 (1954) 

[YamOOa] M. Yamagami, K. Matsuyanagi, Nucl. Phys. A672, 123 (2000) 

[YamOla] M. Yamagami, K. Matsuyanagi, M. Matsuo, Nucl. Phys. A693, 579 (2001) 

[Zhu03a] S. Zhu et al, Phys. Rev. Lett. 91, 132501 (2003) 



